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Sound waves passing through arcs in air are photo- 
graphed with a rotating mirror camera, whereby the sound 
velocity in the arc core is determined. From the observed 
sound velocity the arc temperature may be calculated by 
a method previously described. The arcs measured by this 
method, together with the sound velocities and corre- 
sponding temperatures, are: 


Cored carbon 1.85 X 10° cm sec. 5500°K 
NaCl cored C 1.57 4740 
Al cored C 2.14 6160 
W cored C 2.17 6220 
W 6 mm electrodes 2.04 5950 
W 9.5 mm electrodes 2.27 6440 
W-Fe 2.13 6150 
Fe-Fe welding arc 2.08 6020 


The experimental error corresponds to +100° at 5500°K. 
The experimental requirements of the sound source are 


discussed in detail. The gas temperature across the arc 


core is found to be constant within the error of observation. 
For most of the arcs studied the current density and elec- 
tric gradients have also been measured, so that the data 
may be used as a test of the thermal ionization theory. 
The “‘effective ionization potential (I.P.)’’ of the arc gas 
is in all cases above the I.P. of the metallic vapor and 
below that of nitrogen. For the W arcs the effective I.P. 
is 15 volts, midway between the first I.P. of O2 and No. 
For the arcs other than tungsten, this quantity is lowered 
appreciably by electrode vapor. A small partial pressure 
of Na lowers the arc temperature by approximately 1500°. 
The temperature of a welding arc between a coated iron 
welding rod and a steel plate is measured to be 6000°, from 
which the dissociation of the arc atmosphere is found to be 
81 percent. This high degree of dissociation points to the 
importance of. the dissociation-diffusion-recombination 
process for heat transfer in welding arcs and explains the 
recently reported negative results of attempts to arc weld 
in pure argon. 





1. INTRODUCTION 


N optical method of observing the passage 
of sound waves through arcs has been de- 
veloped and is herewith described. The basis of 
the method is the experimental fact that a sound 
“wave” moving through a radiating gas is visible 
by virtue of the increased luminosity in the high 
pressure region. This wave front, which travels 
ata velccity 5 to 8 times greater than the velocity 
of sound at room temperatures, is photographed 
by means of a high speed camera; from these 
photcgraphs the sound velocity is measured and 
the gus temperature calculated. 

This means of determining the sound velocities 
is superior in simplicity and accuracy to the 
method previously described,':? and is limited 

1C, G. Suits, Proc. Nat. Acad. Sci. 21, 48 (1935); 


Physics 6, 190 (1935). 
?H. Poritsky and C. G. Suits, Physics 6, 196 (1935). 
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only in that very high luminosity is required in 
the arc. 


2. EXPERIMENTAL METHOD 


This experimental method is illustrated in Fig. 
1, where the sound produced by a condensed 
spark discharge in air is passed through an are 
and viewed by a rotating mirror camera. 

As a sound source, a spark discharge in air is 
very satisfactory for this purpose. From an ex- 
tensive literature’ it is known that an oscillatory 
or nonoscillatory spark in air produces a single 
expanding pressure pulse which, for photographic 

3 Mach and Gruss, Wiener. Ber. 78, 14 (1878); 98, 277 
(1888); Mach and Salcher, Wied. Ann. 32, 277 (1888); 
Raps, Wied. Ann. 50, 193 (1893); Toepler, Ann. d. Physik 
14, 838 (1904); Foley and Souder, Phys. Rev. 35, 373 
(1912); Foley, Phys. Rev. 16, 449 (1920); J. E. Smith, 


Phys. Rev. 25, 870 (1925); Okubo and Matuyama, Phys. 
Rev. 34, 1474 (1929). 
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end Film 


Fic. 1. Apparatus for photographing sound waves in arcs. 


purposes, is approximately one mm thick. The 
velocity of propagation of the ‘‘wave’’ front is 
initially high but soon attenuates to a constant 
value of approximately 346 m sec.~ (at 25°C) 
within a short distance from the source. The dis- 
tance required for this reduction to normal 
velocity depends upon the discharge rate, con- 
denser energy, spark voltage, and probably other 
variables. This point is discussed in more detail 


below. 


By reference to the diagram of Fig. 1, it may 


be seen that the sound spark is so placed that the 
electrode axis is parallel to the optical axis of the 
viewing system. It is apparently necessary to 
make the observations along an element of the 
“sound cylinder” since we have failed to obtain 
readable photographs for any other orientation. 
The direction of passage of sound through the 
arc is, however, immaterial. The image of the 
sound wave in the arc is formed on the slit by 
a lens. A second lens forms an image of the slit 
on the cylindrical film which is scanned by the 
rotating mirror running at high speed. The film 
holder is so arranged that 10 exposures may be 
taken on each film. Both lenses have an effective 
aperture of approximately F3. A spring operated 
shutter is placed at the slit and is so dimensioned 
that the time of opening of the slit is one revolu- 
tion of the mirror, which is approximately 1/300 
sec. 

The synchronization of the shutter with the 
sound spark is accomplished with contacts on 
the shutter which operate the primary of an 
ignition coil, the secondary of which is connected 
to a three-sphere gap that triggers the sound 
circuit. The mirror is rotated on the shaft of a 
single-phase 180-cycle motor operated at 300 
cycles. This motor runs within a percent of 
synchronism to provide a speed of 18,000 r.p.m. 
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This speed is sufficient to give the required time 
resolution and cannot be greatly exceeded in the 
present optical system without photographically 
underexposing the most intense arcs. The expo- 
sure calculated from the slit width, dimensions, 
and mirror speed is approximately 1/2,000,000 
sec., which is insufficient to record some arcs 
properly, as, for example, the low current copper 
arc and the low current carbon arc. 

The experimental method described is similar 
to that used by Anderson and Smith‘ and others 
in studying exploded wires, Beams® and others 
in studying spark discharges, and is similar in 
principle to the apparatus used by Schonland‘ 
in the study of lightning discharges. 

In the shadow method of sound ‘wave pho- 
tography first reported by Foley,*? an image of 
the sound wave is formed as a ‘‘shadow”’ cast 
by the high pressure region on the plate. In the 
present experiment the gas through which the 
wave passes is luminous and the wave front is 
observed by virtue of the higher concentrations 
of radiating particles in the high pressure region. 
This effect was first found in exploded wire 
photographs by Anderson and Smith. 


3. RESULTS 


Three photographic enlargements (approxi- 
mately ten diameters), from rotating mirror 
camera films showing sound ‘‘waves’’ in arcs, 
appear in Fig. 2. Figs. 2a and b were taken for 
the cored carbon projection arc (case 2 below), 
Fig. 2c for the iron welding arc (case 8 below). 
At point 1 of Fig. 2a the exposure of the arc 
appears; at 2 the sound gap has started to 
discharge; at 3 the sound front enters the arc 
and passes through the core, leaving the slanting 
line, the slope of which is proportional to the 
velocity of sound. The striations between 2 and 
3 are a result of the intense light from the 
oscillatory sound gap, scattered in the optical 
system. The portion of this light scattered by 
the hot arc gas is negligibly small. The wedge 
shaped structure to the sound “wave’’ at 3 is 
a result of the direction of viewing which, it will 
be remembered, is parallel to the axis of the 


* Anderson and Smith, Astrophys. J. 64, 359 (1926). 
> Beams, Phys. Rev. 35, 24 (1930). 
6 Schonland and Collins, Nature 132, 407 (1933). 














sound gap but perpendicular to the axis of the 
arc. We have taken some photographs from a 
direction parallel to both the spark and arc 
axes, and in this case the wave front appears as 
a line rather than a wedge. In the best records, 


’ of which the above are examples, the definition 
/ is such that the uncertainty in reading the slope 
, of the wave front is about three percent, which 


at 5500°K corresponds to an uncertainty of +60° 
; in the temperature. 

4 This is probably the largest single source of 
. error in the method. Other errors, in calibration, 
B motor speed, and the like, increase the error to a 
. total of about five percent or +100° temperature 
d uncertainty. The spread of the measurements is 
¥ in most cases greater than this, especially with 
the high current welding arcs, which at best are 
y erratic and quasi-stable. This spread must there- 
fore be interpreted as a real temperature fluctua- 


. tion; it would be remarkable if the welding arc 
. temperature did not fluctuate since the current 
e 


density and other quantities vary. An indeter- 
minate amount of error is involved in the extra- 
polation of the equilibrium constants of the air 


r 
. 


dissociation for the temperatures above 5000°K 
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The uncertainty in this extrapolation cannot be 
very large, however, since the asymptotic limits 
to the composition of the air are known. 

This method of measuring arc temperatures 
is simple and rapid. A film with ten traces may 
be taken very quickly, and on the average there 
are about three satisfactory records per film, 
each a rather accurate determination of tem- 
perature and temperature gradients, if any. 


4. SouUND REQUIREMENTS 


We have obtained satisfactory photographs of 
wave fronts in arcs for values of capacitance 
between 0.15 and 1.2 mf at 33 kv, providing the 
discharge rate is sufficiently high. Flat-rolled 
tinfoil paper condensers are generally unsatis- 
factory because of the high inductance and 
resistance. In particular, for 1.2 mf and 33 kv, 
a discharge frequency of 100 kc is satisfactory, 
whereas the same capacitance and voltage at 
30 ke have not produced readable wave fronts. 

Near a source of intense sound the velocity 
may be very much greater than the normal value 
of 0.346 10° cm sec.~! (25°C), and the question 
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TABLE I. TABLE II. 
1 2 3 1 2 3 4 
SOUND VELOCITY* Errec-  I.P. of 
SOUND SOUND IN PROJECTION TIVE L.P. CHIEE 
VELOCITY AT VELOCITY AT CARBON ARC, 5 cm ELECTRON OF ARC METALLIC 
CAPACITANCE 5 cm (25°C) 10 em (25°C) FROM SOURCE Temp. DENSITY GAS ELEMENT 
(mf) (cm sec.~!) (cm sec.~!) (cm sec.~!) ARC °K nX1074 Vie | 
1.2 58510"  .52510° ~—-1.81 K 10° (A) Copper 4050 39* 95* 769 
6 440 390 1.84 (B) Cored projection C 5500 1.59 11.7 6.9 
a 402 .365 1.83 (C) NaCl cored C 4740 y 8.7 5.11 
15 .380 .350 — (D) Al cored C 6160 1.09 13.6 5.95 
OUI Sneed a (E) W cored C 6220 24 15.3 8.1 
wii TTT ESrséiéi‘é 59502514 BI 
Averages of 10 experimental values. (G) 9.5 mm W 6440 
(H) W welding are _ 6150 
arises as to the possible effect of this abnormal  ‘!) — Fe welding sane 


velocity on the measured velocities in the arc. 
It may be shown’ that the temperature gradients 
between the arc core and the surrounding air 
would be expected to lead to increased convexity 
of the wave front and decreased energy density, 
and, as a result, a rapid reduction of abnormal 
velocities. This prediction has now been veri- 
fied in experiments in which sound waves of 
normal and abnormal velocities have been passed 
through a highly reproducible and stable arc 
where the velocities are measured at the arc 
temperatures, 

The results are tabulated in Table I where the 
first column gives the capacitance used at a 
potential of 33 kv in the sound circuit. In 
columns 2 and 3 are shown the sound velocities 
in air at 25°C at distances of five cm and ten cm, 
respectively, from the sound source, measured 
by a cathode-ray oscillograph method.! In the 
fourth column is the sound velocity of the cases 
listed in column 2 as measured in a carbon arc 
under constant conditions. It may be seen that 
the velocity measured in the arc is constant to 
within two percent for sound waves which in 
air at room temperatures have velocities ranging 
between 15 and 70 percent above normal. We 
conclude from this fact that no appreciable error 
arises from the abnormal velocity effect in the 
observed range. 

The velocity of sound method as described 
has been used to measure the temperatures of a 
number of arcs in air at atmospheric pressure; 
a condensed list of these results is given in Table 
II. In all of the cases observed the temperatures 
lie in the range 4000°-7000°K. Because of the 
variety of experimental conditions, the arcs listed 
in the above table will be described in detail 
separately. 





* Included here as a correction to values previously published in 2, 


Where possible, the current density J and 
electric gradient E in the arc column also have 
been measured for the experimental conditions 
used in the temperature determinations. These 
values are also given below. 

It would be desirable to compare the measure- 
ments of arc temperature by the velocity of 
sound method with similar determinations by 
other methods. The copper arc temperatures 
(case (A)) were taken for experimental condi- 
tions nearly identical with those used by von 
Engel and Steenbeck’?’ and may be compared 
with their value of 4600°+350°. The hope of 
taking measurements by the present method of 
identical cases studied by Ornstein® and his co- 
workers by the band spectra method has been 
fulfilled in only one case since most of the arcs 
of sufficient luminosity for sound wave photo- 
graphs have the bands masked by continuous 
spectrum. However, in the case of the Al cored 
carbon arc, the experimental conditions which 
we have used are sufficiently near those reported 
by Ornstein and Brinkman to make the results 
comparable. The temperature in the violet core 
of this are found by the band spectrum method 
is 6000°K+300°K for an arc length of 1.5 cm 
and a current of 7 amperes. The value which 
we have found for substantially the same case is 
6160°K with a spread of +200°K. From these 
two cases it appears that the results obtained 


7 von Engel and Steenbeck, Wiss. Veroff. a. d. Siemens- 
Konzern 10, 155 (1931); 12, 74 (1933); 12, 89 (1933). 

‘See Ornstein, //is Life and Work, for references given 
therein. 

* Ornstein and Brinkman, Proc. Amst. Akad. Sci. 34, 
498 (1931). 
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ARC TEMPERATURES 


by the velocity of sound method are in good 
agreement with temperature measurements ob- 
tained by two entirely different methods. 


(A) Copper arc 4050°K 


The measurements were made on the arc be- 
tween water-cooled copper electrodes in a stabi- 
lizing air stream tube. The results! are: For arc 
lengths between 3 and 12 cm at a current of 6 
amperes, and for a constant arc length of 9 cm 
and currents varying between 3 and 26 amperes, 
the sound velocity is found to be 1.3810° cm 
sec.—'. The (positive column) current density is 
10+2 amperes cm~ throughout the current 
range, while the electric gradient in the arc 
column is 16.5 volts cm at 6 amperes. 


(B) Cored carbon arc 5500°K 


This arc uses a cored carbon cathode® and 
anode” burning freely in air, anode above, at an 
arc length of 1.1 cm (average), and a current of 
45 amperes. The current density is 50 amperes 
cm~ and the gradient 10 volts cm~. 36-mm 
motion picture projectors commonly use this arc 
as a light source. The sound velocity and tem- 
perature, measured midway between the elec- 
trodes (as in all cases described below), are 
respectively 1.85 10° cm sec.—! and 5500°K. 


(C) NaCl cored carbon 4740°K 


In this arc a cored carbon cathode” was used 
as the top electrode, while the lower electrode, 
the anode, was drilled and filled with NaCl. The 
portion of the arc column near the cathode is 
white and distinct from that portion near the 
anode, which has the characteristic Na yellow. 
The observations were made in the Na column. 
At an arc length of 5.6 cm, and a current of 43 
amperes, for which the current density is 85 am- 
peres cm and the gradient 6 volts cm™, we 
find a sound velocity of 1.57X10° cm sec.~, 
corresponding to the temperature of 4740°K. 


(D) Al cored carbon 6160°K 


The lower electrode was 8 mm carbon, drilled 
and filled with Al powder; the upper electrode, 
the cathode, was an 8 mm cored carbon.’ At 43 


* National Carbon Co., 5/16” Orotip cored. 
° National Carbon Co., High-Low Projector. 
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amperes and 4.2 cm arc length, for which the 
current density is 38 amperes cm~? and the 
gradient 7.9 volts cm, we observe a sound 
velocity of 2.14 10° cm sec.—', corresponding to 
an arc temperature of 6160°K. 


(E) W cored carbon 6220°K 


In this arc the lower electrode, the anode, 
was an 8 mm carbon drilled and filled with 
tungsten powder. The cathode was cored carbon.® 
At 40 amperes and 4.2 cm arc length, with a 
current density of 18 amperes cm~ and gradient 
of 17 volts cm™, we find a sound velocity of 


2.17X10° sec! and an arc temperature of 
6220°K. 


(F) Tungsten 5950°K 


Solid tungsten rods of 6 mm diameter burning 
vertically in air (anode above) were used in this 
experiment. With this diameter of electrode it is 
found that 5.5 amperes is the maximum current 
which can be used without developing instability 
in the arc. At 1.6 cm arc length the current 
density is 20 amperes cm~ and the gradient 
22 volts cm. We find a sound velocity of 
2.04 < 10° cm sec.~, corresponding to a tempera- 
ture of 5950°K. 


(G) Tungsten 6440°K 


With 9.5 mm solid tungsten rods as electrodes 
(cathode above), it is possible to use a current 
of 25 amperes with sufficient stability to allow 
rotating mirror photographs to be made. This 
arc is, however, relatively unsteady and but 
few of the exposures were readable. A large 
spread in the measurements, much greater than 
the experimental error, indicates an actual tem- 
perature fluctuation which is characteristic of 
the instability. Fluctuations in light intensity 
and arc voltage" probably correlate with these 
temperature variations. The average sound ve- 
locity is 2.2710° cm sec.~, nearly 7 times the 
value at 0°C. The temperature corresponding 
to the average velocity is 6440°K; the limits of 
temperature fluctuation are 5900°K and 7600°K 
(estimated), the latter being the highest observed 
in this experiment. 


11C, G. Suits, Physics 5, 380 (1934). 
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5. TEMPERATURES OF WELDING ARCS 


These arcs were used under conditions as 
nearly as possible like those employed in prac- 
tical welding operations. The current was sup- 
plied by a conventional welding generator, and 
the arc length and current were adjusted under 
the direction of a welder to values appropriate 
for the electrodes, polarity, and work. These arcs, 
like all welding arcs, are relatively unstable. 
The arc voltage, light intensity, position, current: 
density, and probably other quantities vary in 
an irregular manner. It was necessary to take a 
great many records to obtain suitable data, but 
some excellent exposures were secured for both 
cases (see Fig. 2c). The arc temperature partakes 
in the fluctuation, so that the spread of the 
measurements was great 


(H) Tungsten-iron 6150°K 


For heating, but not depositing metal, the 
welder uses a 3.2 mm solid tungsten rod as 
cathode. At an arc length of 0.5 cm and a current 
of 80 amperes, the current density varies between 
80 and 110 amperes cm~*. The average sound 
velocity is found to be 2.1310° cm sec. and 
the average temperature 6/50°K. The tempera- 
ture fluctuates between the limits 5000°K and 
7000°K (estimated). 


(I) Iron-iron 6020°K 


In this measurement a 3/16” type W-20 coated 
welding rod was used as cathode above a large 
steel plate. The current was 125 amperes at an 
arc length of 1 cm. The current density fluctuates 
between 94 and 405 amperes cm~*. The observed 
average sound velocity is 2.0810° cm sec.~, 
corresponding to an arc temperature of 6020°K. 
The temperature fluctuation lies between the 
limits 4400°K and 7000°K (estimated). 


6. DiIscUSSION OF RESULTS 


For all of the cases here reported the sound 
velocities across the core are sensibly constant, 
as shown by the constant slope on the photo- 
graphic records, which means that the tempera- 
ture across the core of the arc must be constant 
within approximately 120°, which is the uncer- 
tainty in reading the slope on the film. 

The cores of the atmospheric pressure arcs 
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studied are therefore regions of substantially 
constant high gas temperatures. The gradient 
between the core temperature and room tem- 
perature must occur in the photographically 
nonluminous boundary layers of the arc column, 

The data of Table I allow a test of the thermal 
ionization theory of high pressure arcs first pro- 
posed by Compton.” For that purpose the elec- 
tron density has been calculated from 


n=I/eKE, 


where J is the observed current density, e is the 
electronic charge, E is the electric gradient in the 
arc column, and K is the mobility. The mobility 
K is calculated from Wahlin’s experimental 
values of the electron mobility, assumed to vary 
as (T)' and corrected for the effect of dissociation 
on the mean free path. This latter correction is 
numerically three times at a temperature of 
6000°K. The correction for field is negligibly 
small. 

For this electron density 1, the effective ioniza- 
tion potential V;, is calculated as follows: Using 
the Saha equation" it is assumed that 1 is pro- 
duced by a homogeneous gas of ionization poten- 
tial V;. at the measured temperature and at a 
pressure of one atmosphere. 

This effective ionization potential V;, (column 
3, Table II) may then be compared to the 
ionization potential of N2=16.7 volts, of N=14.5 
volts, of O2=14.1 volts, of O=13.6 volts, and 
finally to V,, (column 4, Table II), the ionization 
potential of the chief metallic element in the arc. 
It may be seen that V;, is, within the spread in 
the measurements, less than 16.7 and greater 
than V,, for all cases tested. For the W arcs 
Vi; 15 volts, which is about midway between 
the ionization potential of O2 and Ne. From this 
we may conclude that the metallic vapor from 
the electrodes does not contribute appreciably 
to the ionization in the tungsten arc. The agree- 
ment between V;, and the ionization potentials 
of oxygen and nitrogen is about as good as can 
be expected for the relatively complicated experi- 
mental conditions, and provides strong evidence 
for the essential correctness of theory. 


12K. T. Compton, Phys. Rev. 21, 266 (1923). 

13 Wahlin, Phys. Rev. 27, 169 (1924). 

4 Compton and Langmuir, Rev. Mod. Phys. 2, 123 
(1930). 















In the cases of the other arcs in Table II the 
effective ionization potential of the arc gas is 
decreased by the presence of the low ionization 
potential metallic vapor from the electrodes. 
Judging by the comparative values of V;. and 
V», the effect is not simply related to the ioniza- 
tion potential of the metal. In general it is 
probably necessary to take account also of the 
electrode volatility and chemical properties. 

The observed constancy of the current density 
in the copper arc is compatible with a variable 
gradient on two grounds: (a) the electron mo- 
bility has a dependence on field strength; and 
(b) a temperature variation with current may 
be occurring within the limits of experimental 
accuracy. It is not possible to decide between 
these two cases without greater accuracy in the 
temperature measurements, but it is true that 
the dependence of mobility on field strength is 
sufficient to allow a constant current density at 
constant temperature (and hence m) within the 
limits of measurement of J. | 

The interesting effects produced in the arc by 
the introduction of sodium vapor are: (1) a re- 
duction in arc temperature; (2) an increase in 
current density; and (3) a decrease in gradient. 
The removal of metallic vapor by the use of W 
electrodes produces the reverse effect. This be- 
havior is to be expected on theoretical grounds. 
The observation that the W arc in air is an 
air’ arc raises some questions concerning the 
causes of its high luminosity. The partial pres- 
sure of W vapor is very small, yet the luminosity 
is greater than any other case studied. It is also 
true that the electrodes are rapidly consumed in 
air, and that the spectrum is predominantly 
continuous. In pure hydrogen, however, W elec- 
trodes evaporate very slowly and the discharge 
is practically colorless. 

These facts indicate that the high luminosity 
is due to incandescent tungsten oxide. It is 
probably true that these particles are in equi- 
librium with the hot arc gas, but that is specula- 
tion at present. 


“é 


Welding arcs 


The measurements on the welding arc fix the 
average temperature at 6000° and the limits of 
temperature fluctuation at 4400° and 7000° for 
this type of discharge. The type W-20 coated 
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welding rod arc represents about the most stable 
type of welding arc in common use. 

As already pointed out, a small amount of 
metallic vapor of low ionizing potential lowers 
the arc temperature, increases the current den- 
sity, and decreases the gradient. The first and 
third of these factors contribute to arc stability. 
It is probable that one important function of the 
coating material of welding rods is to provide 
this low ionizing potential vapor. 

From a knowledge of the welding arc tempera- 
ture and recently determined'® equilibrium con- 
stants for the dissociation processes 


O.=2 O 
N22 N 


we may specify the degree of dissociation in 
welding arcs. In Fig. 3 are shown the total partial 
pressures of the atomic fractions in dissociated 
air at atmospheric pressure. At the temperature 
of the welding arc, 6000°K, 81 percent of the 
arc atmosphere is atomic oxygen and nitrogen. 
The effectiveness of dissociated gas tn heating 
metallic surfaces, which act as catalysts for the 
recombination process, is well known'® and be- 


1 Giauque and Clayton, J. Am. Chem. Soc. 55, 4887 
(1933); Johnston and Walker, J. Am. Chem. Soc. 55, 179 
(1933); 55, 191 (1933). 

16 Wood, Proc. Roy. Soc. Al02, 1 (1922); Phil. Mag. 44, 
538 (1922); Langmuir, Gen. Elec. Rev. 29, 153 (1926). 
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lieved to be of importance in atomic hydrogen 
welding. The high degree of dissociation in weld- 
ing arcs in air, however, means that these arcs 
are properly regarded as atomic oxygen-nitrogen 
arcs which depend to some extent upon the same 
recombination process for heat transfer. The heat 
transferred by the diffusion of dissociated gas is 
proportional to f the fraction of the molecules 
dissociated, D the diffusion coefficient, and V 
the dissociation potential. In comparison to 
hydrogen, nitrogen has a smaller coefficient of 
diffusion but a larger heat of recombination. 
The relative magnitude of the heat /y transferred 
by the diffusion of H through He and hy for N 
through N» may be calculated from the ratio: 


hy kfuDu Va. AnCuVu 


hy kfsDnVxy .8\nCuVw 





where & is a constant, \ is the mean free path, 
C is the ‘‘molecular velocity,’’ and the subscripts 
refer to the respective gases. If we assume the 
temperature to be the same for both cases, using 


Vu= 4.4 volts, 
Vy =11.7 volts, 
we find that 

hy/hy =3.3. 
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Thus, when the relative degrees of dissociation 
(fa assumed to be 1), diffusion coefficients, and 
dissociation potentials are taken into account, 
the diffusion heat transfer by atomic nitrogen is 
about one-third as effective as by atomic hydro- 
gen at 6000°K arc temperature. The temperature 
of the atomic hydrogen arc with 1.5 mm electrode 
separation, in the current range 4-8.5 amperes, 
has been measured by Ornstein, Eymers, and 
Wouda,” who find T=4900°—6300°K, with a 
dependence upon current. As shown by Lang- 
muir,!® the heat transfer by direct conduction 
is very small compared to that transferred by 
diffusion, so that we must conclude that the 
diffusion process is of major importance to the 
functioning of the welding arc in air. 
Experimental evidence supporting these con- 
clusions is found in the negative results of recent 
attempts'® to arc-weld in an atmosphere of pure 
argon. It is found that at the same or greater arc 
wattages which produce satisfactory welds in 
air there is a marked lack of pool formation and 
fusion of the deposited metal, indicating an in- 
effective heat transfer. The absence of a dissocia- 
tion process in the monomolecular argon is 
probably accountable for these results. 


17 Ornstein, Eymers, and Wouda, Proc. Amst. Akad. Sci. 
34, 505 (1931). 
18 Doan and Schulte, Phys. Rev. 47, 783 (1935). 
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The shot effect of secondary electron currents from 
plates of triode tubes was measured at 10° cycles. Calcula- 
tions of the electronic charge from shot effects of temper- 
ature-limited primary currents, to check the method, 
agreed with Millikan’s value within 1 percent. Measure- 
ments were made on five tubes having different secondary 
emissivities. The average number of secondaries per 
primary electron, w, ranged from / to 4.5. The variation 
of w with energy of the primaries is shown to agree quali- 
tatively with existing theory. In every case the shot effects 
of both grid and plate currents were measured. Campbell’s 
general expressions for the shot effect of secondary emission 


plate currents were extended to grid currents. From the 
observed differences between grid and plate shot effects it 
is concluded that secondary emission and primary impact 
are simultaneous within 10~® second or less. The large shot 
effects found at high values of w are correlated with the 
fluctuations of w from one impact to another. A qualitative 
analysis of secondary emission processes shows that, from 
these data, information can be obtained as to the number 
of reflected primary electrons and the number of primaries 
which emit no secondaries. The results obtained are in 
reasonable agreement with data obtained by other methods. 





INTRODUCTION 


N their paper on the shot effect, in 1925, Hull 

and Williams! described a single measurement 
on the shot effect of secondary electron emission 
from the plate of a vacuum tube. This measure- 
ment was made at 750 kilocycles under experi- 
mental conditions such that the average current 
of secondary electrons was just equal to the 
average primary current, so that the net plate 
current was zero. They concluded from their 
results that emission of secondary electrons was 
“independent” of the primary impact. The 
primary current in their tube was partially 
space-charge limited. The present paper is a 
report of measurements on a number of tubes 
under conditions in which the complication due 
to space charge was not involved, and with 
secondary currents up to 4.5 times the primary 
current. 


THEORY OF SECONDARY EMISSION SHOT EFFECT 


Stimulated by the work of Hull and Williams, 
and following his own papers on fluctuation 
phenomena in general, Campbell” * has derived 
an expression for the shot effect of a plate current 
consisting of primary and secondary electrons. 


* This work was done 1926-28 in the Research Labora- 
tory of the General Electric Company. See L. J. Hayner 
and A. W. Hull, Phys. Rev. 33, 281A (1929). 
eas) W. Hull and N. H. Williams, Phys. Rev. 25, 147 

*\N. R. Campbell, Proc. Camb. Phil. Soc. 15, 117, 310, 
513 (1909-1910). 

§N. R. Campbell, Phil. Mag. 50, 81 (1925). 


He has considered two cases: (1) that in which 
the emission of secondary electrons is “‘simul- 
taneous” with the impact of the primary, and 
(2) the case in which these events are “‘inde- 
pendent.’”’ The exact significance of the terms 
simultaneous and independent involves the fun- 
damentals of the theory of the shot effect. Fry‘ 
and others have shown that the properties of the 
measuring circuit must be considered in any 
attempt to connect measurements on shot effect 
with the fluctuations of the electron stream. 
They begin with the fact that a single electron 
striking the plate is equivalent to a charge 
placed on the condenser of the plate circuit. 
They then assume that this charge is dissipated 
by an aperiodic or oscillatory current given by 
the ordinary theory of transient currents. Some 
experimental support for this assumption lies in 
the fact that the work functions of metals can, 
in order of magnitude at least, be calculated 
from the classical image force. They then com- 
pute the desired result, the power dissipated in 
the circuit, by summing the effects due to the 
individual electrons. 

_ Campbell considers the effect on the circuit of 
two separate events, (1) the arrival of a primary 
electron, and (2) the departure of a secondary. 
If the transient current produced by the primary 
is completely neutralized by the opposite current 
due to the secondary, the two events are “‘simul- 
taneous.”” The net effect on the circuit is then 
zero. If there is no neutralization, the two 


4T. C. Fry, J. Frank. Inst. 199, 203 (1925). 
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events are “independent.’’ The net effect is 
then the sum of the separate effects. Evidently 
the two events will always be registered as 
simultaneous if the time interval between them 
is small compared with the period of the circuit, 
if oscillatory, or to its time constant, if aperiodic. 
On the other hand, if the interval between the 
events is long compared to the period or time 
constant, then in most cases the events will be 
registered as independent. In one possible but 
improbable case, if the time interval between 
the two events is exactly equal to a small multiple 
of the circuit period, the effect on the circuit 
will be nearly the same as if the effects were 
simultaneous. 

In case there is more than one secondary 
electron, an additional time factor enters, namely 
the time between the emission of the several 
secondary electrons. In Campbell’s derivation, he 
has tacitly assumed that all the several sec- 
ondaries come off simultaneously even when 
they are independent of the primary. The result 
for the case when the several secondaries come 
off independently of each other as well as of the 
primary is also given below. 

It is assumed that the arrival of the primary 
electrons is completely random. There is every 
reason to believe that temperature limited 
currents fulfill this condition. If N is the average 
number arriving at the plate per second, and X? 
the mean square deviation from N, we have, by 
the definition of random, 

X?= N=iv/e, (1) 
where ip is the average plate current and e the 
electronic charge. 

Let each primary electron produce an average 
number w of secondaries. Let the mean square 
deviation from w be 7”. The average number of 
electrons gained by the plate per second will 
then be N(1—w). Let the mean square deviation 
from this number be X,?. Campbell® shows that 
if the emission of the secondaries is simultaneous 
with the arrival of the primaries, the result is 
X,°= (io/e){(1—w)?+7°} (simultaneous). 


(2) 


If the emission of all secondaries is simultaneous, 
but is independent of the time of impact of the 
primaries, then 
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X 2 = (io/e) {1+2+7?} (half independent). 


In addition there is the case not considered by 
Campbell in which the secondaries are inde- 
pendent of the primary and of each other. 


(3) 


X,?=(io/e){1+w} (independent). (4) 


7 has no meaning in this case, since the sec- 
ondary current, like the primary, is completely 
random. 

These equations may be applied directly to a 
system composed of a filament emitting elec- 
trons, surrounded by a grid and a plate. The grid 
is made sufficiently positive to produce saturation 
of the electron stream. The plate is kept at a 
positive potential lower than that of the grid, 
so that secondary electrons emitted by the plate 
are collected by the grid. Secondary electrons 
emitted by the grid do not escape from the grid 
and are not observed. 

Let the average number of electrons emitted 
by the filament per second be N=i)/e. Some 
fraction ¢ of these will be caught by the grid. 
This fraction will depend upon the geometry of 
the tube and the relative values of grid and 
plate potentials. Then (1—¢)N per second will 
strike the plate so that (2), (3), and (4) now 
become 





X P= (1—$)(to/e){(1—w)*+77}, (5) 
X P= (1—$)(io/e) {1 +0? +n7}, (6) 
X = (1— $)(to/e) {1+}. (7) 


The mean square deviation X,? of the current 
to the grid is also of interest. The N@ primary 
electrons collected by the grid are completely 
random and independent of the (1—¢)N pri- 
maries giving rise to secondaries at the plate. 
Hence 

X F= (io/e){o+(1—¢)(w?+ 77) } (8) 


for the simultaneous and half independent cases, 
and 


X = (io/e){o+(1— 4) a} (9) 


for the fully independent case. 

So far, only the fluctuations of the electron 
stream have been discussed. These fluctuations 
ultimately cause a deflection of a meter in the 
output circuit of an amplifier. The value of this 
deflection, of course, involves the properties of 
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the shot effect tube and the circuits associated 
with the tube and amplifier. Several methods 
have been devised for taking these properties 
into account, by means of a comparison between 
the deflection due to the shot effect and the de- 
flection produced by a measured sinusoidal 
voltage introduced into the system at some 
point. The theory of any such method, when 
applied to a pure temperature-limited current, 
gives J, the average electron current, as a func- 
tion of e, the circuit constants, and the measured 
value of the comparison voltage. The comparison 
between theory and experiment consists in com- 
paring the current thus computed with the 
current zp read on a d.c. meter. If similar methods 
are applied to currents which are not purely 
random, or do not involve single electrons, the 
calculated value of J will represent that random 
current whose fluctuations would equal the 
fluctuations actually occurring. The value of J 
thus computed may be less than Zo, as in space- 
charge limited currents, or greater, as in the 
flicker effect. 

To connect the computed values of J with the 
expressions for the fluctuations, it is only neces- 
sary to use Eq. (1) in combination with one of 
the Eqs. (2) to (9). For instance, from (1) and 
(5), and (1) and (8), using J, and J, to dis- 
tinguish the plate and grid shot effects, re- 
spectively, 


J,/e=X?= (1— 9) (to/e) {(1—w)? +77}, (Sa) 
J,/e=X P= (io/e){o+(1—9$)(w*+n7)}, (8a) 


and similar equations from (6), (7) and (9). 

Certain very important conclusions may be 
drawn from Eqs. (5a) to (9a) without any 
knowledge of w or 7’. If the difference J,—J, is 
formed, 


J,—Jp=to{1—2(1—9)(1—w)} 
simultaneous (10) 


Jo—Jn=tolo—(1—¢)} 
half independent or independent. (11) 


The values of found under the conditions of 
this work range from 0.20 to 0.88. The right side 
of (11) may therefore be either positive or 
negative. The right side of (10) is positive if, 
as in these experiments, w>1, and reduces to 


io — 2i,. Here 7, is the average plate current, and 
is to be taken positive if the plate gains electrons; 
its value is read from the d.c. plate meter. 
Hence Eq. (10) may be tested without knowing 
w, 7? or even ¢. It may be stated here that the 
results to be described below confirm (10) and 
exclude (11) as a possibility. 

From (8a) and (9a) 7? may be computed. It 
should be pointed out that in so doing the 
correctness of Campbell’s theory is assumed in 
order to derive information about 7’. 


ng = Jp/(1—$)io— (1—w)? (12) 
from measurements on the plate current. 
10° = (Jo— div)/(1—$)io— w* (13) 


from measurements on the grid current. The two 
values of 7? thus obtained from a given set of 
data should agree. The determination of ¢, w 
and 7? will be discussed below. 


EXPERIMENTAL METHOD 


At the beginning of the work measurements 
were made with a resistance coupled amplifier 
and a shot circuit tuned to 96.6 kc. The method 
of Hull and Williams was followed ; this involves 
introducing the comparison voltage directly 
across the first stage of the amplifier, together 
with a measurement of the equivalent series 
resistance of the shot circuit. The amplifier, 
however, had a linear and nearly flat frequency 
response curve, which eliminated the laborious 
calculations of the amplification factor F in the 
Hull and Williams method. 

The values of e calculated from temperature 
limited shot effects observed with this apparatus 
were always from eight to fifteen percent too low. 
The amplifier was therefore discarded in favor 
of the tuned-circuit amplifier described below. 
Subsequently it was discovered that, due to 
penetration, a correction of about 5 percent to 
the calculated inductance of the inductance 
potentiometer was necessary. This was in such a 
direction as to bring the e values within +2 
percent and —5 percent of the correct value. 
The remaining discrepancies were probably due 
to insufficient shielding. 

A number of shot effect measurements on 


secondary emission were made with this amplifier 
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Fig. 1 Circuit diagram of the shot effect apparatus, showing the shielding in broken lines. A, 
first amplifier tube; S, shot tube; P, inductance potentiometer; A, attenuator; Gi, input galva- 


nometer. 


by using a comparison method. Measurements of 
temperature limited shot effect were made fre- 
quently during the secondary emission observa- 
tions, and the e calculated from the former was 
used in calculating the J’s for the latter. In the 
light of the subsequent discovery of the correc- 
tion factor, the results obtained in this way are 
believed to have about the same precision as 
the later results, and they are therefore include 
in this report. ; 

Most of the shot effect measurements were 
made with a tuned shot circuit of natural fre- 
quency 99.5 ke and a_tuned-circuit-coupled 
amplifier. The sinusoidal calibration voltage was 
introduced into the shot circuit according to 
the method of Williams and Vincent? in which an 
inductance potentiometer of very low impedance 
is connected permanently in series with the 
inductance or capacity of the shot circuit. The 
shot effect equation in this case is 


to = (C*wo"/2e)(v;:"/A), (14) 


where ip is the average temperature limited 
electron current; C is the total capacity of the 
shot circuit; wo/2m the resonant frequency of 
the shot circuit and amplifier; e the electronic 
charge; v,; the sinusoidal voltage which, when 
introduced into the shot circuit, gives the same 
deflection in the output system of the amplifier 
as the mean square shot voltage due to 79; and A 
is an amplification factor defined below. If the 
calibration voltage is furnished by a length / of 
an inductance potentiometer of inductance L, 
per unit length, carrying a current kJ,, where k 


5 N. H. Williams and H. B. Vincent, Phys. Rev. 28, 1250 
(1926). 


is an attenuation factor, then 
io = (CPwoth*IL 1? /2e)(I,?/A). 


Here A=fo°y(f)df, where ¥(f) is the amplifica- 
tion at frequency f, ¥(fo) being set equal to 1. 
¥(f) is obtained from the power output of the 
amplifier when the frequency is varied, keeping 
the current kJ, through the inductance po- 
tentiometer constant. 

The method of Williams and Vincent has 
several advantages. First, any regeneration 
present produces the same effect in the amplifica- 
tion of the mean square shot voltage as in that 
of the calibration voltage. Second, no measure- 
ment of resistance is necessary. Another ad- 
vantage is that the amplification factor A is 
easily calculated. The method has several dis- 
advantages, all of which are connected with 
shielding problems. These will be discussed 
below. 


(15) 


THE APPARATUS 


The amplifier was similar to that described by 
Hull and Williams. There were five stages, the 
first four having 222 tubes and the last a 201A. 
In the last stage high plate voltage and large 
negative grid bias were used to insure linear 
amplification of very large signals. The stages 
were coupled with tuned circuit impedances. 
Suitable chokes and bypass condensers were 
inserted in all battery and galvanometer leads to 
prevent regeneration. 

The output detector was a 700-ohm thermo- 
couple inductively coupled to a coil in the plate 
circuit of the last tube. The thermocouple for 
measuring the input currents (see Fig. 1) was 

















calibrated with d.c. The deflections of its 
galvanometer were accurately proportional to 
the square of the calibrating current. The de- 
flections of the output galvanometer produced 
by the sinusoidal calibrating voltage were exactly 
proportional to the deflections of the input 
galvanometer ; showing that the output thermo- 
couple was a square detector. This is essential to 
the method. The input and output galva- 
nometers were mounted together on a Julius 
suspension. Both zeros were constant within 
0.1 mm. 

The shot circuit and the method of applying 
the calibrating voltage are shown in Fig. 1. 
The capacity C includes the distributed capacity 
of the coil, the capacity to ground of the grid or 
plate of the shot tube, and of the grid of the first 
amplifier tube, and the capacity of all the con- 
nections to ground. Its value was 2318 micro- 
microfarads, +0.1 percent. 

The sinusoidal voltage was supplied by a 
pliodynatron oscillator’ with a frequency of 
99.5 kc. The frequency was measured with a 
General Radio precision wave meter with an 
estimated precision of 0.1 percent. This fre- 
quency was at all times constant within the 
limits of the precision of measurement. 

The inductance potentiometer was similar to 
that described by Hull and Williams and others. 
It was made entirely of brass and very carefully 
constructed. A correction’ of about 5 percent for 
penetration was necessary in calculating Li. 
The value of ZL; was 3.746X10-* henry +0.1 
percent. The same oscillator and inductance 
potentiometer were used with both amplifiers. 

When the calibration voltage is applied in 
series with the shot circuit as shown in Fig. 1, 
the voltage required from the inductance po- 
tentiometer is very small. In this case it was less 
than 1/100 of the voltage which would have 
been required if the calibrating voltage had been 
applied directly on the grid of the first amplifier 
tube. In order that the input galvanometer read- 
ing should not be too small for a desirable pre- 
cision, the attenuator shown in Fig. 1 was used. 
This attenuator, of a ladder type, was con- 
structed of very small high resistance wire. Its 


® A.W. Hull, Proc. I. R. E. 6, 5 (1918). 
7 A. Russell, Phil. Mag. 17, 524 (1909). 
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constant, calculated from d.c. voltage and resist- 
ance measurements, was k= 0.00764+0.00001. 

The amplifier and shot circuit were in metal 
boxes inside a double walled wire cage. The 
oscillator and various parts of the input system 
were in separate metal boxes as shown in Fig. 1. 
Special shielding precautions were necessary 
since this type of input circuit is inherently 
sensitive to stray signals from the oscillator. 

The tuned shot circuit was given additional 
protection by shortening the connecting wires 
and enclosing the inductance in a grounded 
copper box within the usual metal shot box. 

Because of the presence of the attenuator, only 
a small fraction of the measured current through 
the input thermocouple eventually reaches the 
inductance potentiometer. Therefore any cur- 
rents induced in the attenuator will produce 
relatively large errors in the potentiometer 
currents. Since the attenuator is of high im- 
pedance, these induced currents are inevitable 
unless it is well shielded from the coil delivering 
power from the oscillator. With the unshielded 
attenuator the input-output curve was not a 
straight line but a complicated curve not passing 
through the origin. The necessary shielding was 
accomplished by surrounding the attenuator 
with a brass cylinder closed at the attenuator 
end and extending for several inches over the 
end of the inductance potentiometer. This brass 
cylinder was insulated from the potentiometer 
and inner cage, and was connected to the outer 
cage. 

A third source of stray signals is currents in 
the shielding metal near the zero end of the in- 
ductance potentiometer. These cause a signal on 
the amplifier when the oscillator is turned on 
with zero length of the potentiometer in use. 
Also, if signals are delivered to the amplifier 
varying / and keeping J; constant, the output 
will not be proportional to 7. This error is 
especially serious since the undesired voltage 
bears a definite phase relation to the calibration 
voltage, so that correction may not be made by a 
subtraction of the initial deflection from the total 
amplifier output, as in the case of two inde- 
pendent voltages. With the squared output scale, 
an initial deflection of 2 mm in a total deflection 
of 200 mm may mean an error of 20 percent in 
the mean squared calibration voltage. It was 
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found that a stray signal of 4 10-® volt would 
cause an error of 1 percent in the calculation of 
io from (14) when ipo was about 0.5 ma. This 
stray signal was sufficiently eliminated by 
making all the grounds of the amplifier and shot 
circuit, and of all the shielding inside the inner 
cage to the exact zero end of the inductance 
potentiometer. The tests of the effectiveness of 
this system of grounds were (1) that with /=0 
no output deflection was observed with an 
amplification several times that necessary for 
measuring the shot effect; and (2) that for 
different values of / with J, constant, the out- 
puts were proportional to / within less than 1/2 
percent. 


THE EXPERIMENTAL PROCEDURE 


During a series of measurements the effective 
plate resistance of the shot tube varies; this 
alters the equivalent series resistance R of the 
shot circuit. Although R does not appear ex- 
plicitly in Eq. (14), an increase in R decreases 
the amplification at wo, and broadens the reso- 
nance curve, thus increasing both v;? and A. 
The change in amplification at wo is auto- 
matically measured whenever the calibrating 
voltage is applied ; but to redetermine A for each 
experimental point is laborious and time con- 
suming. It was desirable to make the measure- 
ments as quickly as possible since there were 
gradual changes in the sensitivity of the surfaces 
emitting secondary electrons. Hence A was pre- 
determined once and for all as a function of the 
amplification at wo over the whole range of 
resistance change introduced by the shot tubes. 
This was done by inserting small series re- 
sistances, 2.36 ohms to 16.6 ohms, in the resonant 
circuit. After the calibration it was only necessary 
to keep the filament currents and plate voltages 
of the amplifier constant with the aid of meters. 
A direct check of A from time to time showed 
the validity of this procedure. A was found to 
vary less than 0.2 percent with time. 

The half-width of the amplification curve 
ranged from 1.5 to 2.8 kc, and the corresponding 
values of A from 1.75 to 2.98 ke. 

The actual shot measurements were made in 
the following way. Conditions were chosen so 
that the shot effect to be measured gave an 
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output deflection d,; equal to approximately 1/2 
full scale. Then a signal was introduced into the 
shot circuit from the inductance potentiometer 
such that the output deflection was approxi- 
mately doubled, its value becoming dz. Because 
of the random nature of the shot effect d; and d, 
were unsteady. A satisfactory procedure was to 
make ten readings of each at regular time 
intervals. For convenience, instead of attempting 
to make d exactly equal to 2d,, the input current 
was slightly corrected to get the value of J, for 
use in Eq. (15). The filament circuit of the shot 
tube was then opened, and an output d; due to 
the calibrating voltage alone was obtained. The 
ratio d3/(d2—d,) gave the relative amplification 
necessary to obtain A from the curve previously 
prepared. It was necessary to leave a resistance 
in the shot circuit while measuring shot effects of 
secondary emission plate currents because of 
the negative resistance characteristic of the shot 
tube. 

A test of the reliability of the apparatus and 
experimental procedure was made by evaluating 
e from the data on several temperature limited 
shot currents ranging from 0.2 ma to 1.0 ma. 
The values of e ranged from 1.57010- to 
1.604 10-'® coulomb, with an average of 1.587 
<10-'* coul. Calculation of the precision of e 
from Eq. (15), assuming very conservative values 
for the precisions of the various measured 
quantities gave a precision of +1.5 percent. 
In the secondary emission measurements there 
was an added uncertainty due to the slight 
unsteadiness of the secondary emission. It is 
estimated that the precision of the J’s was not 
worse than +2 percent. 


THE SECONDARY EMISSION TUBES 


The shot effect measurements on secondary 
emission were made on five three-element tubes 
designated A, B, C, D, E. The measurements on 
Tubes D and E were made with the resistance 
coupled amplifier. Tubes A, B and C were of 
type 201A, D and E of type 210. These tubes 
were sensitized for secondary emission by Mr. 
H. C. Thompson, except for Tube C, which was 
an ordinary factory tube. The material used for 
sensitization was BaO. The grid wires were of 
nickel, and the plate material was molybdenum 
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Fic. 2. Hyatt’s values of (1—@) as a function of E,/E, for 
a 201A tube. E, =360 volts. 


or nickel. Typical plate voltage-plate current 
curves for similar tubes are shown by Hull.® 
In the tubes used secondary emission began at 
plate voltages between 5 and 10 volts. The plate 
current became zero (w=1) at voltages ranging 
from 20 to 124 volts. These values varied some- 
what with the grid voltage. This is probably due 
to variations of angle of incidence of the primaries 
due to the geometry of the tube. At higher plate 
voltages the plate current became negative, and 
in some cases was over three times the total 
emission from the filament. 

Although it is possible to tell from the plate 
current-plate voltage curves that at the high 
plate voltages several secondary electrons are 
emitted for each primary incident on the plate, 
it is impossible to calculate w unless the fraction 
(1—@) of the primary electrons which strike the 
plate is known. Since (1— ¢) is a function of the 
relative grid and plate voltages, it must be 
measured experimentally. These measurements 
were kindly made for the author on a tube of 
type 201A by Dr. J. M. Hyatt, using the 
published® method of calibration with positive 
caesium ions. The curve showing the relation 
between (1—¢) and E,/E, for E,=360 volts is 
shown in Fig. 2. This curve is decidedly different 
from the curves Hyatt obtained with a tube 


8]. M. Hyatt, Phys. Rev. 32, 922 (1928); 33, 1100A 
(1929). 
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Fic. 3. Relation between yoy and log E, for tubes A, B, 
and C. 


having cylindrical elements and an axial fila- 
ment. In that case, for E,/E, from 0.1 to 2.0 
the fraction of positive ions caught by the grid 
was found to agree closely with the ratio of 
solid to open portions of the grid. In the 201A 
type tube with closed flattened plate and grid and 
hairpin filament, the plate receives in general 
fewer electrons than would be predicted from the 
geometry. It should be noted, however, that 
the fraction striking the plate as E,/E, ap- 
proaches 1 is 0.88, which agrees with the value 
0.87 calculated from the dimensions of the grid 
wires and supports. 

Fig. 3 shows the relation between log w and 
log E, for tubes A, B and C with E,=368 volts. 
Lower values of E, gave curves similar in shape 
but somewhat displaced to the right, corre- 
sponding to the slightly higher values of E, for 
w=0. The curves are all characterized by a linear 
portion at the smaller values of E,, and then a 
region where w increases less rapidly, or even 
decreases as E, increases. The maximum comes 
at higher voltages for the more sensitive tubes. 
All these tubes have the same geometry. It is 
certain that the falling off of w at higher voltages 
is not due to space charge or to errors in @. 

The following equations hold for the linear 
portions of the curves for A, B, and C, re- 
spectively : 


wa = (E,/29)-*; we =(E,/54)-*; wo = (E,/105) 4. 
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Data obtained by Hyatt and Smith® and by 
Copeland'® on heat-treated molybdenum when 
plotted in this way show similar curves. The 
straight portions yield exponents of approxi- 
mately 0.34 and 0.44. Frohlich" has derived 
relations between w and E, for the ends of the 
range of E,. He finds that w should increase in 
proportion to E,'/? at low voltages, whereas 
when £E, is a few hundred volts, w should de- 
crease slowly with increasing E,. 


RESULTS 


Typical values of the secondary emission shot 
effects J, and J, are shown plotted against w in 
Fig. 4. The units of the J’s are, as stated above, 
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Fic. 4. Plate and grid shot effects as a function of w, for 
tube A. E, =368 volts. in =0.303 ma. 


milliamperes of temperature limited shot current 
which would give the same fluctuations as the 
secondary emission currents. The curves shown 
are for tube A with E,=368 volts. The curves 
for different grid voltages are similar, but not 
identical since ¢ is involved in the magnitudes of 
J, and J,. For tube B, which is a poorer sec- 
ondary emitter, the curves are higher, while for 
tube E, the best emitter, the J values are lowest. 
Tubes C and D, for which measurements were 
made only at w=1, also follow this rule. 

The first object of the secondary emission 
measurements was to decide between Eqs. (10) 


* J. M. Hyatt and H. A. Smith, Phys. Rev. 32, 929 (1928). 
10 P, L. Copeland, J. Frank. Inst. 215, 593 (1933). 
tH, Fréhlich, Ann. d. Physik 13, 229 (1932). 
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Fic. 5. J,—J, asa function of i9—2i,, for all tubes and 
voltages. The 45-degree line is the theoretical relation for 
the simultaneous case. 


and (11). It was found that J, was always larger 
than J,, and hence Eq. (11) cannot hold. To 
verify Eq. (10) J,—J, has been plotted for all 
the observed experimeatal points against ip — 27, 
Fig. 5. A 45° straight line representing Eq. (10) 
has also been drawn, and it is evident that most 
of the points lie on this line within the expected 
precision. For instance, the large deviation of the 
largest J,—J, is more than accounted for by 
assuming a 2 percent error for each of the J's. 

In this figure a number of results are invisible 
because the value of 7, was zero, so that all such 
points would fall together on the curve. These 
points are shown separately in Table I. In this 
case Eq. (10) reduces to J,—J,=io. The last 
two columns show the agreement between these 
quantities. The discrepancies are on the average 
about 1 percent. The agreement is thus better 
than would be expected from the precision of the 
J's. This is due partly to the fact that the sec- 
ondary emission was steadiest at the low plate 
voltages. Furthermore, since J, was much less 
than J, the magnification of errors arising from 
taking the difference was less serious. 

The above data show definitely that primary 
impact and secondary emission are simultaneous 
within 10~® second or less, since it is reasonable 


















TABLE I. Summary of secondary emission shot effect 
measurements at ip=0, (w=1), showing the range of grid and 
plate voltage and also the verification of the relation Jg— Jp =1%. 











E E J J Jg—J i 
TUBE (voits) (volts) (ma) (ma) (ma) (ma) 

A 368 28.5 0.357 0.0420 0.315 0.303 
319 30.5 0.349 0.0454 0.304 0.303 
184 32.0 0.381 0.0714 0.310 0.303 
138 33.55 0.395 0.0950 0.300 0.303 
92 35.0 0.440 0.137 0.303 0.303 
B 368 56.0 0.395 0.0978 0.297 0.303 
184 58.0 0.476 0.169 0.307 0.303 
C 368 108 0.503 0.193 0.310 0.303 
230 124 0.548 0.245 0.303 0.303 
D145 48 0.162 0.058 0.104 0.100 
145 48 0.302 0.107 0.195 0.200 
145 48 0.625 0.225 0.400 0.400 
E 230 20 0.230 0.0274 0.203 0.200 
46.0 27.5 0.309 0.110 0.199 0.200 
46.0 275 0.311 0.111 0.200 0.200 








to expect that some deviation from the equation 
expressing simultaneity would appear if the 
average time between impact and secondary 
emission were one-tenth the period of the shot 
circuit. 

The quantity 7? can be calculated from either 
J, or J, by Eqs. (12) and (13). This has been 
done for tubes A, B, C, and E, with Hyatt’s 
values of ¢. These do not apply directly to tube 
E, which is of the 210 type. However, the size 
and spacing of the grid wires are the same in 
the two tubes so that in the range of E,/E, 
where ¢ is nearly constant and equal to the 
geometrical ratio, the ¢’s should be the same. 
The additional precaution was used of keeping 
E,/E, constant at 0.598 for all values of w 
during the measurements on tube E. The values 
of 7,2 and 7,2 plotted against w are shown in 
Fig. 6. 

Like the J’s from which they are computed, 
the values of 7? are larger for the poor secondary 
emitters. At w=1, 7? ranges from about 0.5 to 1 
for the tubes used. The values of 7” in general 
increase more rapidly than w. The fact that 
n, is somewhat different from 7,2 in the lower 
portions of the curves is again due to the 
multiplication of errors arising in the com- 
putation of 7,?. The values of 7,? are considered 
much more reliable. The curves in Fig. 6 include 
measurements made with various grid voltages. 
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Fic. 6. 7? as a function of w, for Tubes A, B, C, and E, 


and for all voltages. The circles refer to plate fluctuations 
and the crosses to grid fluctuations. 


The relation between 7? and E£, is linear within 
the experimental precision for a given tube and 
grid voltage. The slope and intercept of the line 
vary with £, and are also different from tube to 
tube. 


DISCUSSION OF THE RESULTS 


Theoretically the probable time interval be- 
tween primary impact and secondary emission 
can be estimated in several ways. The com- 
parison of computed and measured conductivities 
of metals leads to the conclusion that the mean 
free path of an electron entering a metal is of 
the order of 100 atom diameters.” This means 
that an electron entering a metal with 100 volts 
energy will lose its excess energy and become a 
member of the normal Fermi distribution in a 
time of the order of 10-™ second. The time 
interval between primary and secondary electron 
should be of this order of magnitude also. 

Another estimate comes from the time during 
which the incident electron may be expected to 
stay in an excited level inside the metal (neg- 
lecting collisions). The probability of an electron 
transition to a given lower level in a metal should 
be of the same order as that for a similar transi- 
tion in an atomic system. Because of the enor- 
mous number of available energy states in the 


2 A. Sommerfeld and H. Bethe, Handbuch der Physik, 
second edition, Vol. 24/2, p. 354. 
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metal, lifetimes of the order of magnitude given 
above should result. 

It is therefore not surprising that the results 
indicate simultaneity, and some other explana- 
tion must be found to account for the result of 
Hull and Williams. 

It is of interest to examine what would be 
expected theoretically for 7*. To do this it is 
necessary to consider the known facts of sec- 
ondary emission produced by primary electrons 
of several hundred volts energy." * © 

(1) At all incident energies there are some 
secondary electrons which have velocities nearly 
equal to that of the primaries. These are equiva- 
lent to reflected electrons, and there can be only 
one of them produced by any one primary. 
If such electrons are produced by a fraction x of 
the primaries, the value of w for this fraction is 
equal to 1. Their contribution to the plate shot 
effect is zero. The value of x may be expected to 
depend on the nature of the surface and on the 
speed of the primaries. From the papers of 


Farnsworth" it appears that for moderately 


clean, pure metal surfaces x varies from 0.5 or 
more for 10-volt primaries down to a few percent 
for 100-volt primaries; and that x may be even 
greater for surfaces having adsorbed gas layers. 

(2) There are some primaries which produce 
no secondaries because they communicate all 
their momentum to ions in the metal (inelastic 
collisions), or to electrons which are directed 
inward, or to electrons whose outward velocity 
component is not sufficient to carry them out 
through the surface. The number of these 
primaries will thus depend on the velocities and 
angular distribution of the secondaries inside 
the metal. Both factors may themselves depend 
on the energy of the primaries. Scattering of 
primaries and secondaries must be taken into 
account. All these considerations make the 
estimate of this fraction very difficult. If it is 
assumed, however, that the metal atoms scatter 
both primaries and secondaries like elastic 
spheres, the chances are relatively small that no 
secondaries would get out. Let the fraction of 
primaries which produce no secondaries be y. 


13 W. Bothe, Handbuch der Physik, first edition, Vol. 24, 
y ? 

r 4H. E. Farnsworth, Phys. Rev. 20, 358 (1922). 
6 T, Soller, Phys. Rev. 36, 1212 (1930). 
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Then w,=0; and the contribution to the plate 
shot effect, J,, is that of an equivalent tempera- 
ture limited current. 

(3) Of the remaining primaries, some will 
produce many secondaries, some few. If this 
fraction is z, then J.=[(1—w.)?+n2]. The 
range of values which w, and 7,2 may assume is 
predictable. Since each secondary loses an energy 
equal to or greater than the work function on 
its way out of the metal, the maximum number 
of secondaries which can be produced by any 
one primary is (E,+a)/(E.+a), where a is the 
work function, £, the primary energy, and EF, 
the average secondary energy, all in volts. If 
a=3 volts for the BaO surfaces used in this 
work, and E,=4 or 5 volts from the data of 
Farnsworth, the maximum number becomes 3 to 
40 for primary energies from 20 to 300 volts. 
There may also be as few as one secondary: per 
primary. The cases in which this happens 
cannot be distinguished from cases of reflection 
except by experiments on velocity distribution 
such as those of Farnsworth" and Soller." 

The quantity 72 may have any value from 
zero to w,’, depending on the distribution of 
the number ejected about the mean value w,. 
It can be shown that neither of these values will 
fit the data for any possible values of x, y, and z. 
Because of the large number of random factors 
involved in the secondary emission process, the 
most likely assumption is that the number 
ejected varies in purely random fashion. In this 
case 7, would be approximately equal to w, 
itself, at least for high primary voltages. 

It is now necessary to derive the relation 
between w, and 72 and the w and 7? of Campbell’s 
less detailed picture. For unit primary current, 


yt2(1—w,)=1-—w 
or 
w.=(at+w—1)/s if (y+s)=a, (16) 
then 


Jp={y+s(1—w.)? +2} =(1—w)?+n°. (17) 
Putting in nz =w, and the value of w, from (16), 
2? =w(1—w) +(a+w—1)?/z. (18) 


It is therefore possible to determine certain 
limits for a and z from the experimentally ob- 
tained values of w and 7’. First of all, when 
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w=1, Eq. (18) becomes a?/z= 7’, or as an upper 
limit a=7*. For tube A this would mean 40 
percent reflected electrons, since 7’ =0.6; for the 
less sensitive tubes B and C, a would be higher. 
The reason why at w=1 there is less reflection for 
the less sensitive tubes is that the primary 
voltage is higher. These results are in agreement 
with the data of Farnsworth and of Soller. 

For large values of w, (a—1) will be small 
compared to w. Eq. (18) then gives 


P=ot+w(1/z—1). 


The data for the highest value of w for each tube 
then lead to 2=0.96+0.01 for tube E; z=0.82 
+0.04 for tube A, and z=0.5+0.2 for the tube B. 
In the latter case, the largest value of w was 
only 2.8, so that the assumption (a—1)<a, is 
not good. In the middle ranges of w little can be 
said of the parameters a and z. 

The fact that at high primary energies z is so 
large for some tubes shows (1) that there is 
very little reflection, and (2) that there must be 
considerable scattering; for otherwise most of 
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the secondaries would receive inward momenta, 
and the fraction y would be large. 

The author takes pleasure in acknowledging 
her great indebtedness to Dr. A. W. Hull, who 
suggested this research and who contributed 
valuable suggestions and advice throughout the 
progress of the work. She also wishes to express 
her appreciation to Dr. J. M. Hyatt for the 
measurements of ¢, and to Professor Bernhard 
Kurrelmeyer for assistance in taking the data 
and for discussion during the preparation of the 
paper. 

Note added in proof August 23, 1935. In recent 
work the author has found that the variation of 
w with grid voltage, observed throughout this 
work, is probably due mainly to secondary elec- 
trons from the grid, a few of which will have 
sufficient energy to reach the plate. Such elec- 
trons may also increase the grid and plate shot 
effects by amounts of the order of ten percent 
at the higher plate voltages. The relation 
J,—J»=%—2i, should remain accurately valid. 
The effects due to grid secondary electrons will 
be discussed in detail in a future paper. 
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Measurements have been made of the equilibrium per- 
meabilities of softwoods to air of different relative vapor 
pressures. Transverse sections less than the average fiber 
length in thickness, in which the open cavities account for 
practically all of the permeability, show practically no 
change in permeability with changes in the equilibrium 
relative vapor pressure. Sections thicker than the maximum 
fiber length, the permeabilities of which are dependent 
upon the size of the pit openings, show an increase in per- 
meability with a decrease in equilibrium relative vapor 
pressure. When the square root of the permeability is 


plotted against the moisture content of the wood in equi- 
librium with the various relative vapor pressures of air, 
practically a linear relationship is obtained from 0 to 20 


_ percent moisture content. At higher moisture contents the 


permeabilities are considerably less than the linear rela- 
tionship calls for. This is due to films of water forming 
across the capillaries. Higher pressures than those used 
are required to overcome the effect of the surface tension 
of the water in these capillaries. A new means of deter- 
mining the distribution of size of openings in a porous 
membrane based on these findings is given. 





HE capillary structure of wood has been 
studied from the pulping standpoint by 
Johnston and Maass® and by Sutherland, John- 
ston, and Maass," and from the preservation 
standpoint by Buckman, Schmitz and Gortner.? 
The author has determined the effective capillary 
sizes in water-saturated wood by a combination 
of pressure permeability and electro-osmotic flow 
measurements.” * Electrical conductivity meas- 
urements”: '° that were also used have since been 
shown to be incapable of giving capillary sizes 
because of the erroneous assumption that the 
conductivity of the cell walls is negligible. 

The only available data to give an indication 
of the effect of shrinking and swelling of wood 
upon the effective capillary sizes are those of 
Buckman, Schmitz and Gortner® which indicate 
that the effective capillaries increase in size upon 
drying of the wood. More complete information 
along this line should be of considerable value in 
determining and understanding seasoning, pres- 
ervation, and pulping practice. 


EXPERIMENTAL PROCEDURE 


The equilibrium permeability of wood sections 
to air of different relative vapor pressures was 
measured with a modification of the differential 
pressure drop apparatus previously described 
(see Fig. 1).8 It consists essentially of three com- 

*Presented before the Colloid Division, American 
Chemical Society, at its 89th meeting in New York City, 
April 22-26, 1935. 


+ Maintained at Madison, Wis., in cooperation with the 
University of Wisconsin. 


partments, 1 and 2 separated by the wood section 
and 2 and 3 separated by the standard capillary 
tube C. Manometers were provided to determine 
the pressure drop between compartments 1 and 2 
and 2 and 3. The means of holding the sections 
was modified as follows: Holes were drilled in one 
end of cylindrical pieces of wood W in the fiber 
direction with a centerless bit, the bottom serving 
as a transverse section. A soft rubber stopper R 
connected the section to the apparatus. It 
showed no tendency to loosen or exert an appre- 
ciable stress on the section over a .complete 
moisture change cycle. The wood section was 
enclosed in a large cylindrical glass tube with 
rubber stoppers at each end. A stopcock S was 
provided in compartment 2 to prevent the 
passage of air through the wood section during 
the course of bringing it to moisture equilibrium. 
This was done to cause humidification to occur 
from the outside of the specimen so as to more 
nearly duplicate normal drying and reabsorption 
conditions. An electric heating coil was wrapped 
about the exit tube B to prevent condensation 
of water from air at a high relative vapor pressure 
when being discharged from the apparatus. The 
apparatus was held at a temperature of 40°C in 
a thermostatically controlled water bath. The 
air was humidified by bubbling through satur- 
ated salt solutions in glass jars filled with glass 
beads. About an inch of mercury was placed in 
the bottom of the jars and the air first bubbled 
through this to prevent salt from working back 
into the jet. The air passed through a trap filled 
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with glass wool to remove any traces of entrained 
spray before entering the apparatus. The relative 
vapor pressures given in Table I were obtained 
from gravimetric determinations of the moisture 
content of the air at approximately the same rate 
of flow as that used in the experiments. 

This apparatus proved satisfactory for deter- 
mining the permeability of transverse sections. 
It could not be used, however, for radial and 
tangential flow as practically all the air would 
pass through the end-grain part of the cylindrical 
wall of the wood section rather than through the 
end. In the case of the transverse sections the 
passage of air through other than the end-grain 
ends was entirely negligible (1 percent or less). 
Further, this glass apparatus could not be used 
for other than the transverse sections because it 
would not stand the pressures required for ob- 
taining measurable permeabilities with some of 
the resistant sections. It was thus necessary to 
devise another means of making the measure- 
ments using a more rugged apparatus. 

It was found that sections clamped in the old 
type of flush rubber-gasketed flange, when care- 
fully marked, could be removed from the clamp 
and reinserted in practically the same position 
with a variation in the pressure drop readings of 
not more than 2 to 3 percent. It was hence 
decided to use this type of flange clamp and 
remove the sections and prehumidify them to 
approximately the equilibrium relative vapor 
pressure of the air to be used in the next perme- 
ability test. In this way checking of the sections 
was entirely avoided. Although the source of 
error in these measurements was somewhat 
greater than in those made with the apparatus 
of Fig. 1, there was the distinct advantage of 
being able to make measurements on one series 


TABLE I. Relative vapor pressures obtained with different 
saturated salt solutions. 


CAPILLARY STRUCTURE OF WOOD 











RELATIVE 
VAPOR 

MATERIAL TEMPERATURE PRESSURE 

°C 

Distilled water 40 96.3 
Saturated BaCle 40 86.7 
Saturated NaCl 40 73.4 
Saturated MnCl2 40 50.5 
Saturated MgCle 40 31.0 
Saturated LiCl 40 11.7 
5 40 0.0 
Distilled water 23-25 98.0 
Saturated BaCle 23-25 89.3 
Saturated NaCl 23-25 74.5 
Saturated MgCle 23-25 31.0 
P20; 23-25 0.0 















































Fic. 1, Low pressure permeability apparatus used for 
making measurements on transverse sections of wood. 


of sections while another series was being pre- 
humidified. 

More substantial humidifiers were also neces- 
sary for these measurements. They were made of 
galvanized iron pipe with fritted glass disks 
clamped to the lower end of the pipe to break up 
the air into fine bubbles. To prevent the salt 
from clogging the fritted glass disks, air was 
always passed through a water prehumidifier 
first. The salt solutions thus always reduced the 
relative vapor pressure of the air rather than 
raised it. These humidifiers were used only at room 
temperature, 23° to 25°C. The relative vapor 
pressures obtained with them are given in Table 
I. The sections were prehumidified in humidity 
rooms held at 80°F and 30, 75, 90 and 97 percent 
relative humidity, respectively, and in a desic- 
cator over P.O; for at least a week. The equi- 
librium relative vapor pressure change in the 
apparatus was so small that readings could be 
made as quickly as pressure equilibrium was 
established. 

The velocity of the flow of air through both 
the standard capillary tube and the wood sections 
can be expressed by Poiseuille’s equation: 


V.=areéP./8yl. 
(for standard capillary tube c), 


Vw = NOQrrtP/8nylw (for the wood section), 


(1) 
(2) 











in which V,. and V,, are the velocities of flow 
through the capillary and through the wood 
section, r. and 7, are the radius of the capillary 
and the average effective radius (to be more 
exact, the fourth root of the average fourth 
power) of the wood capillaries, /. and /,, are the 
corresponding lengths and P. and P,, the corre- 
sponding pressure drops, 7 the viscosity of the 
air, N the number of effective capillaries in the 
wood sections in parallel per unit of cross section 
and Q the effective cross section, all expressed in 
centimeter-gram-second units. When the stand- 
ard capillary and the wood section are connected 
in series as in this investigation 


V.=V. (3) 
and 
Nr.‘/le=réP./Q Pe. (4) 


Both N and r, are unknown so that in order to 
determine either, further data are necessary.* ° 
During the course of the shrinking and swelling 
of wood, N will, however, remain constant unless 
checking occurs. The pressure drop ratio P./P,, 
and thestandard pressure drop ratiorP./Q/-P 
should thus be proportional to the fourth power 
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of the capillary radii. The change in size of the 
capillary openings in wood have hence been 
expressed in terms of the pressure drop ratio in 
this investigation. 


EXPERIMENTAL RESULTS 


Measurements of the pressure drop ratio for 
thin transverse sections less than a fiber length 
in thickness of the heartwood of initially green 
western hemlock and air-dry white pine at dif- 
ferent equilibrium relative humidities are shown 
in Table II. Air was passed around the outside 
of the sections and out through stopcock S’ 
(Fig. 1) for two days before passing it through 
the section in order to bring them to equilibrium 
with the vapor pressure of the air. Pressure drop 
readings were then made over a period of two 
days. The total pressure drop used was never 
greater than 20 cm of water. The data show that 
the fiber cavities change but very slightly in size 
with shrinking and swelling of the wood. The 
mean deviation of the effective capillary radius 
for the white pine was about 0.5 percent and for 
the western hemlock specimens still less. In the 
case of the white pine a part of the permeability 


TABLE II. Effect of changes in relative vapor pressure upon the permeability of transverse sections of wood less than a fiber 
length in thickness. 








DIMENSIONS OF SECTIONS 





SPECIES THICKNESS CROSS SECTION 


CAPILLARY RELATIVE VAPOR EQUILIBRIUM 
CONSTANT PRESSURE PRESSURE MEAN 
reA/le OF AIR DROP RATIO DEVIATION 





Cm Cm 
White pine heart 0.22 1.37 


Western hemlock heart .22 5.30 


Western hemlock heart .40 5.30 


Percent 

1.85 X 1078 0.963 
.867 
.734 
.505 
.310 
-117 
.000 
117 
.310 
.734 
.963 
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CAPILLARY STRUCTURE OF WOOD 


may have been through resin ducts. These are 
entirely absent from western hemlock so that the 
permeability of these sections was entirely 
through open fiber cavities. This approximate 
constancy of the size of the fiber cavities is in 
agreement with findings of Schwalbe and Beiser® 
and Beiser,! and with the density-shrinkage 
relationships: for small sections of wood dried 
under as nearly stress-free conditions as possible 
developed by the author." 

Similar measurements were made on initially 
green transverse sections of both the sapwood 
and the heartwood of western hemlock greater 
than the maximum fiber length in thickness. The 
results obtained for several relative vapor 
pressure change cycles are shown in Fig. 2. The 
curves are numbered in order and the desorption 
part of the cycle indicated with open circles and 
the adsorption part with black circles. There is a 
definite increase in the pressure drop ratio with 
decreasing equilibrium-relative vapor pressure 
indicating that the effective capillaries in the pit 
membranes increase in size when the wood 
shrinks. During the first cycle true equilibrium 
was not obtained or a slight checking of the 
sections occurred. After this the relative vapor 
pressure cycles were quite reproducible. Adsorp- 
tion values of the pressure drop ratio were in all 
cases slightly greater than the corresponding 
desorption values. These hysteresis loops are 
entirely comparable with those for the moisture 
content-relative vapor pressure” and electrical 
conductivity-relative vapor pressure relation- 
ships..* The electrical conductivity-moisture 
content relationship calculated from the rela- 
tionships of each with the relative vapor pressure 
gives a relationship free from hysteresis effects. 
The same is true when the pressure drop ratio- 
relative vapor pressure relationship is converted 
to a pressure drop ratio-mojsture content rela- 
tionship basis. The pressure drop ratio is propor- 
tional to the fourth power of the effective radius 
(see Eq. (4)). The square root of the pressure 
drop ratio is hence proportional to the effective 
cross-sectional area of the capillaries. The rela- 
tionship between the square root of the pressure 
drop ratio and the equilibrium moisture content 
is plotted in Fig. 3 for the same data as in Fig. 2. 
A linear relationship free from hysteresis effects 
is obtained below a moisture content of approxi- 
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Fic. 2. Effect of changes in the relative vapor pressure of 
air passing through transverse sections of the heartwood 
and the sapwood of western hemlock upon the resulting 
equilibrium pressure drop ratio. Numbers indicate order 
of measurements: Open circles, desorption; closed circles, 
adsorption. 


mately 20 percent. The external volumetric and 
cross-sectional swelling of wood is directly pro- 
portional to the moisture content from oven-dry 
to the fiber-saturation point." The cross-sectional 
areas of the pit membrane openings thus appear 
to be inversely proportional to the extent of 
swelling of the membranes. Evidently the thin 
pit membranes respond to changes in the relative 
vapor pressure more rapidly than does the rela- 
tively heavy rim about the membrane. In the 
course of desorption the membrane tends to 
shrink on drying. As the response of the heavy 
rim is considerably slower the membrane is put 
under tension. This tension is relieved by internal 
shrinkage which results in an increase in the 
cross section of the openings. 

Similar determinations of the relationship 
between the moisture content and the square root 
of the equilibrium pressure drop ratio were made 
on sections of both seasoned white pine heart- 
wood and initially green heartwood and sapwood 
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Fic. 3. Effect of changes in the moisture content of 
transverse sections of the heartwood and the sapwood of 
western hemlock upon the square root of the equilibrium 
pressure drop ratio. Numbers on curves refer to those of 
Fig. 2. Open circles, desorption; closed circles, adsorption. 


of western hemlock cut in the three different 
structural directions using the second apparatus. 

The measurements were made with total 
pressure drops ranging from 5 to 100 cm of 
mercury for the different sections. Readings at 
at least three different pressures were taken with 
each section. The pressure drop ratio P./ P. was 
plotted against the total pressure drop P. plus 
P., and the value of P./P, extrapolated to zero 
applied pressure was used in the calculations. 
This graphically corrects for impact turbulence 
effects.’ The pressure drop ratios built up to a 
maximum and then decreased to an equilibrium 
value which was used in these calculations. This 
effect for the vapor flow was considerably less 
than that reported by Buckman, Schmitz and 
Gortner? for liquid flow as might be expected for 
electrokinetic effects would undoubtedly be less 
for vapor flow than for liquid flow. 

The results of these measurements with the 
second apparatus are given in Figs. 4, 5, and 6. 
In all cases the square root of the pressure drop 
ratio increases in a linear manner with a decrease 
in the moisture content from about 20 percent 
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Fic. 4. Effect of changes in the moisture content of 
transverse sections of the heartwood of white pine upon the 
square root of the equilibrium pressure drop ratio. Open 
circles, desorption; closed circles, adsorption. 


to oven dry. Further, the ratio of the square roots 
of the pressure drop ratios at 0 percent and 20 
percent moisture content are practically constant 
for sections of different thickness and even for 
sections cut in the different structural directions, 
but the values seem to vary with the species. 
The values of this ratio are given in Table III 
together with the thicknesses of the sections and 
the standard pressure drop ratio for unit effective 
thickness. The measurements made by the first 
method are for a varying cross section caused by 
shrinking and swelling of the section. The meas- 
urements made by the second method are for a 
fixed cross section determined by the size of the 
opening in the flange. As the wood shrinks the 
number of capillaries that are effective for flow 
increases. For this reason the ratios of the square 
roots of the pressure drop ratios obtained by the 
second method are higher than those obtained 
by the first. The values obtained by the second 
method were corrected to the basis of the first 
by dividing by 1 plus the cross-sectional shrink- 
age from 20 to 0 percent moisture content (see 
last column of Table III). 

These ratios of the square roots of the equi- 
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CAPILLARY STRUCTURE OF WOOD 


librium pressure drop ratios seem to depend . 
upon the nature of the membranes traversed 
rather than the arrangement and number og ° 
traversed in series or parallel. If the pit mem- 
branes which are composed of lignin have a sk 
similar hygroscopicity to that of the wood as a < 
whole, which seems to be the case, the volumetric 
shrinkage of the membrane substance from 20 to 
0 percent moisture content should be about 20 
percent. If it is assumed that the rim of the 
membrane changes dimensions by a negligible 
amount during the course of the shrinkage of the 
membrane, hypothetical values for the ratio of 
the square roots of the pressure drop ratios for 
0 and 20 percent moisture content can be cal- 
culated for various specific cases. If the shrinking 
manifests itself entirely as a change in the 
thickness of the membrane, the ratio would be 
1.12. If the shrinkage of the membrane manifests ‘Pr | nes! 
itself entirely by increasing the size of the open- 4 \> 
ings in the plane of the membrane, the value of ie Osc ~~ 

the ratio would depend upon the fractional void 
cross section of the membrane. If this were one- 
half, the ratio would be 1.20, if one-quarter, the inorstone ConrenT Gencenry” 
ratio would be 1.40. In reality there is very likely af 0.5, Pet of changes in the moisture content of 
a change in both the thickness of the membrane pine upon the square root of the equilibrium pressure drop 
and in the size of the openings. If one-half of the ratio. Open circles, desorption; closed circles, adsorption. 
shrinkage manifests itself in the direction of the ratio of the square roots of the pressure drop 
thickness of the membrane and the other half ratios will be 1.26. This hypothetical value is in 
in the plane of the membrane and the void cross reasonable agreement with the actual experi- 
section of the membrane is one-quarter, then the mental values. 
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TABLE III. Changes in permeability of wood with changes in thickness and moisture content. 








RATIO OF THE 
SQUARE ROOTS OF 





THE PRESSURE CoLuMN (6) 
STANDARD PRESSURE DROP RATIOS FOR CORRECTED 
DROP RATIO FOR DRY WOOD AND FOR 
KIND ‘OF SECTION THICKNESS EFFECTIVE UNIT THICK- WOOD AT 20% EXTERNAL 
SPECIES SECTION NO. OF SECTIONS NESS OF DRY WOOD MOISTURE CONTENT SWELLING 
(1) (2) (3) (4) (5) (6) 
Cm 

Western hemlock Sap transverse 10 1.300 8.15 «1078 1.175 1.175 
Heart transverse 10 1.310 5.25 X 1078 1.200 1.200 

1 919 1.26 X 1078 1.282 1.176 

2 1.857 1.18 X 1078 1.302 1.193 

3 2.880 1.141078 1.288 1.180 

Sap tangential 1 .184 5.57 X 10718 1.215 1.153 

2 .189 1.79 X 10718 1.230 1.168 

Sap radial 1 -186 9.50 X 10712 1.233 1.197 

Heart tangential 1 -196 4.58 X 10718 1.245 1.177 

Heart radial 1 169 6.29 K 10718 1.222 1.183 

White pine Heart transverse 1 1.148 4.13 X 10-9 1.420 1.322 
2 2.125 3.78 X 1079 1.475 1.373 

3 2.980 3.85 X 1079 1.475 1.373 

Heart tangential 1 .158 2.45 X 10712 1.400 1.330 

2 -202 1.56 X 10712 1.385 1.315 

3 .255 1.03 X 10712 1.368 1.300 

Heart radial 1 -112 1.10 x 10-13 1.365 1.335 

2 -140 5.00 x 10714 1.352 1.322 








1 These measurements were made with the first apparatus; r-4/l, was 3.73 X10~* and the effective cross section 1.37 cm’. 
The other measurements were made with the second apparatus; rc//- was 5.27 X 10710 and the effective cross section 0.353 cm}. 
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Fic. 6. Effect of changes in the moisture content of transverse, tangential, and 
radial sections of both the heartwood and sapwood of western hemlock upon the 
square root of the equilibrium pressure drop ratio. All desorption curves. Open 
circles, ordinates as given; half-circles, 10 ordinates; closed circles, 100 


ordinates. 


If the capillary structure of wood were quite 
uniform the standard pressure drop ratios per 
unit of effective thickness should be fairly 
constant for each structural direction. This is 
only approximately true in the fiber direction 
and not at all true in the other directions (Table 
III). All the measurements made on transverse 
sections were for sections only a few fiber lengths 
in thickness. The effective thickness of these 
sections is hence appreciably less than the actual 
thickness as air passing through the sections 
traverses on the average one-quarter of an open 
fiber length at each end. The number of pit 
membranes traversed in series is hence approx- 
imately proportional to the thickness of the 
section minus one-half of the average fiber 
length (0.17 cm). This effective thickness was 
used in the calculations of the standard pressure 
drop ratio per unit effective thickness for the 
transverse sections. The number of membranes 
traversed in series for the other sections was so 
great that the difference between actual and 
effective thickness was negligible. The standard 
pressure drop ratios per unit of effective thick- 
ness decrease somewhat with increasing thickness 
even for the transverse sections. This is due to 
the fact that the probability of maximum-sized 
openings occurring in series decreases with an 


increase in the number of membranes traversed 
in series. It can be readily demonstrated that the 
permeability is greater if the maximum-sized 
openings occur in series rather than when they 
occur in parallel. 

The data show an appreciable deviation from 
the linear moisture content-square root of the 
pressure drop ratio relationship at moisture 
contents exceeding 20 percent. The square roots 
of the pressure drop ratios are less than the 
linear relationship calls for. This is due to the 
fact that moisture condenses in part of the effec- 
tive capillaries even below the equilibrium satura- 
tion pressure because of their minute size and 
eliminates them as a source of flow. The rela- 
tionship between the curvature of a drop and 
its equilibrium vapor pressure is given by Kelvin’s 
equation 


r=2¢M/pRT In po/p, (5) 


in which r is the radius of the drop or of the 
capillary in which it forms, fo» is the saturation 
vapor pressure at the absolute temperature 7, 
p is the vapor pressure over the drop or liquid 
in the capillary, R is the gas constant, @ is the 
surface tension of the liquid, 1 the molecular 
weight, and p the density, all expressed in cen- 
timeter-gram-seconds. The relative vapor pres- 
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7. Distribution of the size of the capillaries effective in controlling the flow of air ina 


tangential section of white pine (7). 


sure in equilibrium with the moisture content at 
which the deviation from the linear relationship 
begins is practically 0.9 for all of the curves. 
This corresponds to a capillary radius of 10 mu. 
It is slightly less than the average effective 
capillary radius previously determined by the 
electro-osmotic method.* At moisture contents 
in equilibrium with higher and higher relative 
vapor pressures more and more of the capillary 
structure is obstructed due to condensation of 
films of water that cannot be broken at the 
pressures under which the measurements were 
made. Measurements were also made of the 
pressure required to overcome the effect of 
surface tension of water in the saturated sections 
using the second apparatus. This prdéved to be 
more sensitive to the detection of the initial flow 
of air than the method previously used. For all 
the transverse sections, detectable displacement 
of water by air occurred at lower pressures than 
previously reported.* These pressures in the case 
of the transverse white pine sections were less 
than the maximum operating pressures used in 
making the measurements on the same sections 
at lower moisture contents (15 to 20 cm of 


mercury). Finite values of the square root of 
the pressure drop ratio were hence obtained at 
the saturation moisture content of 29 percent 
(Fig. 4). In the case of the transverse sections of 
western hemlock (Fig. 6) the pressures required 
to overcome the effect of surface tension ranged 
from 0.5 to 2 kg per square centimeter. Although 
these values exceed the pressures used in the 
measurements, the corresponding capillary sizes, 
70 to 2800 my, are so large that the vapor 
pressure depression by them is less than 0.1 
percent. The square root of the pressure drop 
ratio will thus be zero at practically the satura- 
tion moisture content (30.7 percent). In the case 
of the tangential and radial sections of both 
species the pressure required to overcome the 
effect of surface tension ranged from 20 to 25 kg 
per square centimeter, corresponding to capillary 
sizes of 75 to 60 muy, and relative vapor pressures 
of 98.6 to 98.2 percent. These correspond to 
moisture contents of 27.2 percent for the white 
pine and 28.9 percent for the western hemlock. 

The distribution of size of the capillaries 
effective in controlling the rate of flow of air can 
be calculated from the deviations from the linear 
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moisture content-square root of the pressure drop 
ratio relationship. The deviation of the pressure 
drop ratio or velocity of flow from the theoretical 
value at different moisture contents in equi- 
librium with various relative vapor pressures was 
converted to a pressure drop change-capillary 
radius basis with the use of Eq. (5) for a tan- 
gential section of white pine. The values are 
plotted in Fig. 7, together with the distribution 
curve which was obtained by graphically deter- 
mining the slope of the former curve. Sufficient 
experimental values for the pressure drop ratio 
at different equilibrium relative vapor pressures 
were not obtained to determine the distribution 


curves with any certainty but the example given 
gives the general nature of the curve. The most 
probable radius and average radius are practi- 
cally the same and have a value of about 28 my, 
This value is but slightly larger than the average 
values previously obtained by combining electro- 
osmotic data and pressure permeability data.® 

This method for obtaining the distribution of 
size of pores should be applicable for all kinds of 
membranes in which the openings are sufficiently 
small to appreciably reduce the relative vapor 
pressure. In the case of nonswelling membranes 
the linear part of the relationship would, of 
course, be parallel to the abscissa axis. 
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In previous papers we studied the effects of forces, 
produced by metabolism, on the mechanical stability of 
the metabolizing system. In this paper various other 
effects of those forces are discussed. First of all it is shown 
that, whenever those forces are of such a nature as to make 
the system unstable and cause its eventual division, their 
effect is also that of increasing the size of the system. This 
reminds of the relation between rate of growth and division 
in living cells. Further the effects of those forces on 
molecules and colloidal particles, present in the system, 
but not directly participating in metabolism, are studied. 
It is shown that those forces produce a nonhomogeneous 
distribution of those molecules and particles, and thus 
alter the structure and the physical constants of the system. 
For single molecules the effect is negligible; but for col- 
loidal particles of a size above 10-* cm it may be very large 
and may produce a concentration of those particles at the 
surface of the system, altering thus the permeability of 
the latter. Those effects are present only as long as the 
system metabolizes. This may have a bearing on the 


sudden change of permeability of living cells after death. 
Since, under the influence of those forces the permeability 
of a system becomes a function of the intensity of the 
metabolizing processes, the fundamental equations de- 
scribing such systems are much more complicated than 
those studied hitherto. One particularly interesting feature 
is that those equations now possess in general not one, but 
several stable solutions, so that the specification of all the 
parameters, which determine the external conditions, does 
not determine the configuration of the system in a unique 
way. As shown previously, such a situation leads to various 
hysteresis phenomena. One type of such a hysteresis, par- 
ticularly interesting from the point of view of possible 
biological applications, is illustrated on a numerical ex- 
ample. It consists in an irreversible transition of the 
system from a state of slower growth and greater stability 
into a state of more rapid growth and lesser stability, 
under the influence of certain reversible changes in the 
environment of the system. 





N a series of previous papers,'-* we have 

shown that every metabolizing system is a 
seat of various forces, partly of osmotic, partly 
of atomic origin. Those forces are proportional 
to concentration gradients, produced by metab- 
olism, and are therefore present only as long 
as the latter persists. We have investigated the 
effect of those forces on the mechanical stability 
of the system, and shown how they can cause a 
spontaneous division of the latter. In the present 
paper we shall study some other possible effects of 
those forces. 

An important general effect of the forces 
produced by metabolism is their influence on the 
variation of the size of the system. As we have 
seen, substances diffusing inward produce os- 
motic forces, also directed inward, and tending to 
compress the system, while substances diffusing 
outward produce forces tending to expand the 
system. It is these expanding forces that produce 
under certain conditions, previously studied, a 
spontaneous division of the system. But these 


1N. Rashevsky, Physics 5, 374 (1934). 
* N. Rashevsky, Physics 6, 35 (1935). 
*N. Rashevsky, Physics 6, 33 (1935). 
*N. Rashevsky, Physics 6, 117 (1935). 


same expanding forces will tend also to increase 
the size of the system. Such an increase cannot 
go on indefinitely without some additional pro- 
duction of the substances, of which the system is 
composed. The expanding forces cannot therefore 
be the only factor, which would produce increase 
of size. But it is one of many possible other 
factors, which acts in the sense of increasing the 
size of the system. Other conditions being equal, 
the stronger the expanding forces, the more 
rapidly the system will increase in size. On the 
other hand, as we have seen,‘ the greater the 
expanding forces, the stronger will be the tend- 
ency of the system to divide spontaneously. In 
the particular case, studied previously, of a 
complex metabolizing system, which decomposes 
glucose and partly oxidizes it, we have seen that 
the stronger the glycolysis, the less stable the 
system is mechanically, and the more readily it 
divides. But from what has just been said, the 
stronger the glycolysis, the more rapidly the 
system must also increase in size. This also agrees 
with observations on living cells, the only natural 
metabolizing systems so far known. Rapidly 
growing cells of embryonic tissues are charac- 
terized by a high rate of glycolysis. Still more 
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rapidly growing tumors have a still higher 
glycolytic coefficient. The exact mechanism of 
increase of size of metabolizing systems, de- 
pending on several factors, will be studied 
mathematically in a separate paper. 

Another important effect of these forces is 
their influence on molecules and molecular 
aggregates, present in solution in the system, 
but neither diffusing in or out from the latter. 
Under those conditions their concentration will 
be uniform, unless this uniformity is disturbed by 
some force. 

Let us first consider the ‘“‘atomic forces”’ only. 
The same considerations as apply to the inter- 
action between the molecules B of a nonuni- 
formly distributed dissolved substance and a 
molecule A of the solvent,! hold also for the 
interaction between B and any other molecule C, 
dissolved in the system. In the most general case 
we should investigate also the mutual interaction 
of the molecules B. But this would mean that 
we cannot apply to B the laws of dilute solutions, 
and for the present we shall restrict ourselves only 
to those. However, even for a dilute solution, we 
do in general consider the interaction between 
the molecules of the solvent and the solute.® 

Similarly, retaining the laws of dilute solutions 
for B, we may consider interactions between B 
and C. 

The force, acting on a molecule C placed in a 
field of nonuniform concentrations of B, is given 
by Eq. (8) of a previous paper :! 


4nrnKN 
g= 
M(n—1)(n—3)d"~4 





(1) 


grad c=a grad ¢, 


the notations being the same. If, instead of a 
single molecule, we consider a molecular complex, 
or colloidal particle, consisting of m molecules C, 
then the force will be 


G 4rKNm 
M(n—1)(n—3)d"~‘ 





grad c= ma grad c. (2) 


With® n=6 and d~10-° cm, we find K~10-*, 
This gives: 
a~10745, (3) 
5M. Planck, Treatise on Thermodynamics (London, 
Longmans, Green and Co., 1927). 


*A. Gyemant, Grundzuge der Kolloidphysik (Vieweg, 
Braunschweig, 1925). 
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Consider, for sake of definiteness, that the force 
g is directed from the center to the periphery of 
the system. The work necessary to bring a par- 
ticle from the periphery to the center is equal to 


0 0 
v={ Gdr=af grad cdr=ma(é—c°), (4) 


where é and c® denote the concentrations at the 
surface and at the center, respectively. We have :! 


é—c°=gro/6D. (5) 


The ratio a of concentrations of the particles C 
at the surface and in the center, produced by the 
above forces, is given by Boltzmann’s law: 


a= e~ VED, (6) 


In order that this ratio should be appreciably 
different from 1 we must have: 


y~eT. (7) 


With values for g, 7o and D used previously, we 
find: €—c®°~10-> g-cm™; combining this with 
(3) we find y~mX10-*°. Eq. (7) is therefore 
equivalent to 


m~10°kT ~ 10°. (8) 


(7) and (8) show, that the effect of the forces 
here discussed on single molecules is negligible, 
for then y/kT ~10~° and a is very near to unity. 
For a colloidal particle consisting of 10° mole- 
cules, that is having the size of ~10~* cm, a will 
be already of the order of 3. This means that 
there will be a very strong concentration of the 
particles C near the surface of the system. For 
particles of the order of magnitude of 10-5 cm 
m~10® and a~e!®°, This means that practi- 
cally all particles C will be concentrated at the 
surface, forming a separate surface phase. In 
other words the forces due to concentration 
gradients will under certain conditions produce 
a radical change in the structure of the surface 
of the system. And this change of structure 
naturally must result in a change of permeability 
h. If the system stops metabolizing, or speaking 
biologically, dies, the forces disappear, and the 
particles C, concentrated at the surface, become 
distributed uniformly through the system. This 
is accompanied by a change in the permeability. 
It is known in biology that the death of a cell is 














accompanied by almost instantaneous changes in 
its permeability. 

The osmotic forces, which are of the same 
order of magnitude as the other ones,! will 
produce similar effects. If V is the volume of the 
colloidal particle, then the total force acting on 
it is: 


(RT V/M) grad c, (9) 


which is formally similar to (2). The calculation 
proceeds in the same way, and leads to similar 
results. 

Of course the force may also be directed 
inwards and then produce a concentration of 
colloidal particles in the center of the system. 
This may result in the separation of a phase. In 
a complex system, consisting of many metabo- 
lizing substances and many different colloids, 
some particles will go one way, some the other. 
But in any case, the concentration gradients due 
to metabolism will produce a heterogeneity even 
in those components of the system which may 
not directly take a part in the metabolic 
processes. 

In a previous paper,? we have shown that 
when the metabolized substances are ionized, the 
system will be the seat of electric charges, those 
charges disappearing when metabolism ceases. 
An excess of electric charges of any sign in the 
system has the effect of producing a “diffuse 
double-layer”’ at the surface.* For charges of the 
order of magnitude of those calculated before® 
the drop of potential ¢ across the double layer 
is of the order of magnitude of 10-* volt.* The 
effect of such a double layer upon the permea- 
bility of the surface to ions is readily seen. A 
positive ion approaching the surface from the 
positive side will pass without other resistance 
than that due to the finite coefficient of diffusion. 
If the ion approaches the surface from the 
opposite side, in order that it should pass 
through, we must have (1/2)mv? >3-3 X10 ge, 
m being the mass of the ion, v the normal 
component of its velocity, e the charge expressed 
in e.s.u. and g being expressed in volts. In other 
words the permeability will be asymmetrical. The 


permeability / of a surface is defined by 
n=h(ci—C2), (10) 


where » denotes the number of molecules, 
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passing per unit area per unit time, when c; and 
c2 are numbers of molecules per unit volume on 
each side of the surface. We can therefore cal- 
culate the permeability of an electric double- 
layer for an ion in the ‘‘direction of smaller per- 
meability,’’ by putting c.=0 and calculating the 
fraction n/ci=h of the ions which can pass 
through the layer, that is the fraction for which 
(1/2)mv? >3-310- ge. This fraction is cal- 
culated in a similar way, like in the old theory 
of thermionic effect’ with an electronic gas inside 
the metal, and a double layer at the surface.’ In 
this way we find for the number of ions passing 
through in the “difficult direction’ per cm? per 
second the expression 


ny’ =¢\(k/2am)'T'e—/T (11) 

with b=3-3X10-*ge/k. (12) 
The corresponding permeability is equal 

hy=n,/c1. (13) 


The permeability of the other direction is 
obtained by putting b=0. This permeability he 
is also equal to the permeability in the same 
direction for ions of opposite charge. The ratio 
of the two permeabilities is equal to 


hy/he=e~?'T, (14) 
With g~10- volt we find b~10?. 
hi/he~1/e. (15) 


A slight increase in g causes a very large decrease 
of hi/he. For ¢=50 millivolts, h;/h2=0.04. 

This ‘elective permeability” for ions of one 
sign remains only as long as the electric charges 
remain, that is as long as the system meta- 
bolizes. As a rule living cells show an elective 
permeability for ions, which disappears after 
death. The present technique of permeability 
determination in single cells, which is largely 
qualitative, does not give any evidence for an 
“asymmetry” of the kind here discussed. But 
such an asymmetry for ions is not improbable in 
view of the fact that some biological membranes 
show an asymmetry even for neutral molecules. 

As a general result we find that metabolic 


7E. Bloch, Thermionic Phenomena (Methuen and Co., 
London, 1927). 
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processes affect by a purely physical mechanism 
the permeability of the system for the meta- 
bolized substances. Since this permeability in 
general may affect itself the metabolic processes, 
we arrive at a generalization of the previous 
treatments. We shall next have to consider h 
and D in our fundamental equations not as given 
constants, but as depending themselves on the 
distribution of concentrations which they effect. 

The study of such general cases leads to a 
number of separate problems of various mathe- 
matical complexity. Mathematically each case 
will depend on the assumptions which we make 
about the nature of the colloidal substances, and 
about their influence on the permeability of the 
surface. These assumptions will determine how 
the permeability will vary with the concentration 
gradients, and since the latter in their turn de- 
pend on the rate of production or consumption 
of a substance, we shall thus study various 
possible laws of variation of / with the total rate 
Q of metabolism of a given substance. Amongst 
many possibilities a particularly interesting case 
is that when the h—Q curve is S-shaped or 
inverted-S-shaped ; that is for small Q, / is nearly 
constant, then in the neighborhood of a certain 
value Q* it more or less suddenly increases or 
decreases, approaching, however, again an asymp- 
totic value for very large Q. A curve of that type 
can be characterized by the initial value /, the 
final value h’, the value Q*, at which the curve has 
its inflection point, that is for which |dh/dQ) is 
a maximum, and finally by the “half-width” AQ, 
which gives the distance along the Q-axis between 
the two points on the curve, for which |dh/dQ)| 
has half its maximum value. AQ is a measure of 
the steepness of the curve. 

The detailed study of this case, which can be 
shown to be rather general, as well as of other 
possible cases, requires rather elaborate mathe- 
matical investigations, which must be reserved 
for separate papers. Here we shall confine our- 
selves to an illustration on an example of a par- 
ticularly important characteristic of such cases. 

Let us again consider a system in which 
glucose is partly glycolized, partly oxidized. 
However, instead of assuming, as we approxi- 
mately did before, constant rates of reactions, we 
shall consider, in better agreement with actual 
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facts, that the rates of reactions are proportional 
to the concentrations of the reacting substances, 

For glucose we shall then have the differential 
equations and boundary conditions: 


DiV2c,=ac,; V%e;'=0; 
hy(ey’ —c¢1) = Dy’ (dey'/dr) = D,(de,/dr) 
for r=7o, (16) 


the subscript ; referring to glucose, the super- 

script ’ referring to the external concentration. 

a is a constant characteristic of the reaction. 
For lactic acid we have: 


DWN?c2=a01— fc2c3; Veo’ =0; 
ho(ce’ —C2) = D2! (dc2' /dr) = D.(dc2/dr) 
for r=fo, 


(17) 


f being another constant, and the subscript , 
referring to lactic acid, and ; to oxygen. Finally 
for oxygen we shall have: 


D3V*c3=1-O7fcec3; V2cs=0; 
h3(c3’ —c3) = D3’ (dc3' /dr) = D3(dc3/dr) 


for r=fo. (18) 


We shall have a similar equation for COz., 
Because of the nonlinearity of the system (16), 
(17), (18) its general treatment presents great 
difficulties. We may, however, simplify consider- 
ably our problem and obtain a general semi- 
quantitative survey of the solution by using a 
device employed previously * namely, we may 
consider with a sufficient approximation the 
case that all D’s are, though finite, very large. 
Then the variation of concentrations inside and 
outside the system will be small as compared 
with the absolute values of the concentrations, 
and in calculating those values we may assume 
them to be uniform. We cannot assume the D’s 
to be infinite, because then the gradients would 
vanish, and with them the forces which we are 
studying. But with finite, though sufficiently 
large D’s, the gradients may be sufficient to 
produce the necessary forces, yet be small 
enough to justify the calculation of the concen- 
trations on the assumption of their uniformity. 
A more rigorous justification of this procedure 
requires a detailed mathematical investigation 
of the general problem and shall be given else- 
where. . 

With this assumption the concentrations may 
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be calculated by the consideration that just as 
much substance enters (or leaves) the system as 
is consumed (or produced) in it.*» * Our system 
of Eqs. (16)—(18) then becomes equivalent to 


4aro2hy(cy’ —€1) = (4/3) ereacy, (a) 
Anrethe(co—C2') = (4/3) rre?(aci— fcec3), (b) (19) 
Anrehs3(c3’ —c3) =1-07(4/3) rre® fcecs, (c) 


where ¢,’, C2’ and c;’ are the given constant con- 
centrations of glucose, lactic acid and oxygen in 
the external medium. To (19) must be added the 
equation for CO:: 


Anriha(cy—c4’) =1-47 (4/3) rre® fcocs. (20) 


The coefficients 1-07 in (19c) and 1-47 in (20) 
are determined from the stoichiometrical char- 
acteristics of the reactions considered. 

The first Eq. (19) gives: 


€,=3hyc1' /(aro+3h,). (21) 


Solving the third Eq. (19) with respect to c3, and 
substituting into the second, we find after some 
rearrangements: 


Cy=— (3h2/aro)co+ [ fhscs’ / (ahs +0 ° 36aro fc) ‘ee 
—3hece'/aro. (22) 


If 0-36roftx<hs, (23) 
then (22) reduces to 
qj = [3h2/arot+fes'/a \co—3hece' /aro. (24) 


Taking c; and ¢2 as coordinates in a plane, let 
us plot the curves, represented by (21) and (24). 
If all h’s are constant, then (21) represents a 
straight line, parallel to the cs-axis, while (24) 
is another straight line, with a positive slope, 
intersecting the c-axis on the positive -side. 
Hence the two lines will intersect in one point. 

If (23) is not satisfied, then (22) represents a 
line, which is slightly convex upwards for small 
values of c2, but which with increasing C2 ap- 
proaches asymptotically the straight line (24). 
There is again one and only one point of inter- 
section between (21) and (22) for positive values 


*N. Rashevsky, Physics 1, 143 (1931). 


of c, and ce. The corresponding values of c; and 
C2 represent the unique solution of our problem. 
Substituting those into the third equation (19) 
and into (20) we find also unique values for cs; 
and ¢4. 

Things change, however, when the h’s are 
not constant, but depend on the total rates of 
consumption of any of the substances, which 
rates determine the small but finite concentration 
gradients. Let us assume here, as an illustration 
only and for sake of simplicity, that the h’s 
decrease with increasing Q3, the total consump- 
tion of oxygen, according to an inverted-S-shaped 
curve discussed above, and are not influenced by 
either Q; or Qe. 

We have: 


Q3=4arrhs(c3’ —Cs). (25) 


Similar equations hold for Q: and Qs. 
Eq. (19) gives: 


C3=h3c3' /(0-36rofcte+hs), (26) 


which introduced into (25) expresses Q3 as a 
function of hs and cz: 


Q3=4- 5Sro3 fcs’hsco/(0-36rofceths). (27) 


On the other hand we have the equations: 


hi=fi(Qs); he= f2(Qs) ; hs=f3(Qs), (28) 


which express the permeabilities as a function of 
Q3. Eliminating Q3; from (27) and each of the 
Eqs. (28) we shall obtain three equations con- 
necting ji, he and hs with ce. From those we can 
express /;, he and hz as functions of ce. If we 
introduce those expressions, which will depend 
on the assumption made about f1(Qs3), f2(Qs) 
and f3(Q3), into (21) and (24) we see that now 
those two equations do not represent any more 
straight lines, but more or less complicated 
curves, which may intersect at more than one 
point. 

The analytical treatment of the problem 
would involve, even in the simplest case, rather 
elaborate calculations. A graphical solution is, 
however, very easily obtained. In Fig. 1, 7 
represents the graph of Eq. (21), @ that of Eq. 
(24) for variable h’s, with the following choice of 
constants : 
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ro=3X10%cm; a=0-2; 
ce’ =5 X10 g-cm™; 
hy=3 X10 cm: min.~!; 


he=3-75 X10 cm-min.—; 
Q32* =4x10-" g -min.~! ; 
hs=4X10- cm-min.'; 
Q33*=4X10-” g-min.~; 


The constants are chosen so as to have biologi- 
cally plausible values. Q3:*, Qs2*, Q33* are the 
values of Q;* for the inflection point of hy, he 
and hs correspondingly. AQs3:, AQz32: and AQ33 
are the corresponding half-widths. We see that 
the two curves intersect at three points Ai, Ae 
and A;. For all points of the curve 7 Eq. (19a) 
is satisfied. All points in the ¢,ce-plane, which 
lie above 7, correspond to such values of c; and 
C2 for which the left-hand side of (19a) is smaller 
than the right-hand side. This means that for 
any such points there is less glucose penetrating 
into system than there is consumed. Hence 
dc;/dt<0 and c, will decrease until a point on 
the n-curve is reached. For all points below 7 
the reverse holds true. c; increases until 7 is 
reached. For all points above @ the left-hand 
side of (19b) is smaller than the right-hand side, 
hence dc2/dt>0O and ce increases. The reverse 
holds true for all points below the curve @. 
Those considerations show that the points A, 
and A; represent stable solutions of (19), whereas 
Az is unstable. A slight displacement in the 
¢:¢2-plane from the point A» will result in such 
variations of c; and cz that either A, or A; will 
be reached. 

We find from Fig. 1 and Eq. (25) and (29) 
that for Ai: 
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f=40; 
c3' =10-* g-cm™. 

hy’ =1-6X10- cm-min.—; 

Qu* =4-4X10-” g- min. 
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cy =10* g-cm-; 


AQs1 =0-8x«10-" g -min.~', 





he’ =0-9X 10-4 cm-min.—'; (29) 
AQs2=1-4X10-2 g-min.—, 
h;’ =2X10~ cm-min.—'; 

AQ33=1-4X10-” g-min.—, 
Q;=13X10- g-min.—; 
Q2=9-4X10-” g-min.—'; 
Q;=3-7X10-" g-min.'. B=0-9 

for A3: (30) 
Q,=10-3X10-” g-min.“; 
Q2=5-6X10-" g-min.—'; 
Q3;=5-2X10-2g-min.— g=0-4 


where 6 denotes the glycolytic coefficient. We 
see that in A, the relative production of lactic 
acid is higher and 8 much larger. 

We thus see that whereas in the previously 
studied cases of constant / the physico-chemical 
configuration of the system was determined in a 
unique way by the parameters which character- 
ize the external medium (such as cy’, ce’ and ¢;’), 
now we find that the system possesses in general 
several configurations of equilibria. We have 
already made previously’ a study of the general 
properties of such systems, and seen that they 
are characterized by various hysteresis phe- 
nomena in which a reversible change of the 
external conditions produces an_ irreversible 
change in the system. Let, for instance, the 
system have the configuration A;. If we reduce 
sufficiently c;’ (the supply of oxygen), then 
because of a decreased Q; all h’s will increase to 
the corresponding maximum values, 7 and @ 
will become practically straight lines, inter- 
secting at the point A’. (Fig. 1, broken lines.) 
If we again restore the original value (10~ 
g-cm-*) for c;’, we shall find that the system will 
pass into configuration A;, and not into the 
original A3, because A’ lies in a region of the 
€1¢o-plane for which dc;/dt<0 and dc2/dt<0. 
The same thing will happen if, with the system 
originally in A3, we temporarily reduce or make 
zero f. In configuration A; the system is char- 


®N. Rashevsky, Zeits. f. Physik 53, 102 (1929); 59, 562 
(1930); 61, 511 (1930). 
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acterized by a lower permeability and a lower 
glycolytic coefficient. In A; both are increased. 
The situation represents a remarkable analogy 
to the fact discovered by O. Warburg,’ that 
temporary asphyxiation of living cells may 
cause them to pass into a condition of higher 
glycolytic coefficient, and higher permeability, 
characteristic of tumor cells, and that this con- 
dition remains even after normal oxygen supply 
is reestablished. 

The above example is used only as an illus- 


10Q. Warburg, The Metabolism of Tumors (Constable 
and Co. London, 1930). 


tration. Similar cases of hysteresis are obtained 
also with dh/dQ;>0. Of course such “multiple 
equilibria’ do not always occur in cases of vari- 
able h’s, but their occurrence is rather general. 
We may consider also the more general case that 
the glycolysis is partly reversible, so that part of 
the lactic acid is transformed back into glucose. 
This complicates somewhat the calculations, but 
leaves the general result unchanged. 

In the more general cases, when we consider 
also the role of the diffusion coefficients, still 
more complex conditions, with more than two 
stable equilibria, may occur. All this will be 
made a matter of future studies. 
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